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ABSTRACT 

The present work investigates the aerothermal effects on the 
electro-optical performance of a model infrared seeker missile at 
Mach 4 and small angle of attack. Two window dome configurations are 
studied: an optical hemispherical and ellipsoidal shape window 
attached to a 3:1 ogive nose. A three-dimensional thin-layer 
compressible flow solver is used to compute the flow field in front 
of the window and the window aerodynamic heating. Grid adaptation is 
performed in order to accurately compute the flow field and capture 
the bow shock. Noise-Equivalent-Temperature-Difference degradation 
due to "background" noise of the hot window is derived and computed 
for 3 levels of array non-uniformity for a range of mean dome 
temperatures covering the whole supersonic regime. Acquisition range 
for a staring sensor is estimated for several maritime scenarios 
modeled by LOWTRAN code. Background photon flux distribution on an 
array generated by the hot dome is computed using radiation transfer 
methods and the computed dome temperature. It is found that array 
non-uniformity has a strong influence on the seeker performance. 
This study indicates that the non-uniformity of the dome temperature 
field has a significant influence on the array fixed pattern noise 
induced by dome emission and, as a result, on the seeker detection 
and tracking ability. 
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I. INTRODUCTION 


Flight at high supersonic speeds through the atmosphere is 
characterized by high aerokinetic heating rates and gas compression 
in front of the missile. Large density variation and heating of the 
window significantly affects the performance of a vehicular infrared 
(IR) sensor, leading to major design problems of IR seeking 
missiles. Presently, IR windows that operate in the low supersonic 
flight regime exist,and are in operation. However, future IR and 
dual mode tactical missiles will require IR systems which are 
capable of effective operation in much more severe flight 
environments. The durability and survivability of the dome plays a 
major role in defining the flight envelopes of these missiles. 

For the case of a classical IR seeker installed at the tip of 
the missile nose, the outermost effect of the flow field on the 
optical performance is due to the large density jump across the bow 
shock wave in front of the IR window. This density gradient refracts 
optical rays and acts as a lens. Gas behind the shock is highly 
compressed and heated. Density variations in this compressed flow 
downstream of the shock create refractive index variations and have 
the effect of a gradient index lens. The principal optical effects 
are resolution (blur and defocus) and boresight error. Another 
effect of the compressed flow is aerodynamic heating of the window 
itself. As the IR dome heats up, its self-emission may adversely 
impact the operation of the seeker system in a high-speed flight 
situation. The added background-induced noise will lower the 
effective performance of the system in terms of detection range. 




resolution and tracking capability. High resolution is desirable 
for detection (including clutter rejection), tracking accuracy, and 
recognition. The boresight error can affect performance in several 
ways. Boresight fluctuations (caused by turbulence, shock 
interaction, or coolant mixing) can blur images and induce angle 
noise in the tracking system. Changes in boresight caused by body 
motion (e.g. changing shock position or window distortion) appear as 
position errors or tracking rate errors that adversely affect 
sensors used in guidance and navigation. 

Significant radiation from the compressed gas comes from heating 
by its passage through a strong shock wave. This effect, however, is 
only significant at higher Mach numbers (= M >7) , when the perfect 
gas assumption breaks down and chemical reactions must be 
considered. 

A mission analysis [1] for the future tactical missile 
generation revealed a requirement for high supersonic speeds of 
operation. This will produce very high aerodynamic pressures and 
thermal stresses on the optical window that protects the IR 
detection system. Severe environmental conditions encountered at 
very high aerodynamic speeds has forced the development of new 
exotic materials, capable of withstanding such high temperatures and 
erosion. Traditionally, the shape of most IR windows or domes that 
protect the IR seeker has been a truncated hemisphere. Although this 
configuration has some advantages from the optics point of view, 
because of easier fabrication and design, it imposes a considerable 
aerodynamic drag penalty on the vehicle [2] . Therefore, an 
ellipsoidal window is proposed. Such window design would decrease 



the aerodynamic drag of the vehicle, increasing its range. In 
addition, it is expected that it will delay boundary layer 
transition and reduce aerothermodynamically imposed stresses within 
the IR window material [3]. As a result, the self-emission photon 
radiation from the window into the detector will decrease [2] [4], 
and the detection range will improve. 

The temperature field on the window is obtained numerically, 
using the OVERFLOW flow solver. Comparison with available data was 
conducted to validate to some degree the computation, and to give 
confidence in the temperature results obtained. Using this data, the 
aero-thermal effects which affect the imaging performance of an IR 
seeker were evaluated. 
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Figure I-l. Infrared Dome Design Constraints 
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II. BACKGROUND 


A. BLUNT BODY STRUCTURE 

A typical blunt-body, supersonic flow field is shown in Fig.2-1. 
This figure shows the general structure of the flow field in front 
of the window mounted at the nose of a supersonic vehicle. The shock 
wave is detached and normal shock relations may be used to estimate 
the change in density along the stagnation streamline. A complete 
flow field and shock-structure solution depends on body geometry and 
flight Mach number. Such solutions are complex because of mixed 
subsonic (iyf<l) and supersonic (M>1) regions in the flowfield. The 
governing equations are mathematically elliptical in the subsonic 
region and hyperbolic in the supersonic region. Time-dependent 
numerical techniques are used to determine flow fields surrounding 
blunt bodies [5]. An empirical formula for the shock stand-off 
distance. A, at the tip of a conical body with a hemispherical nose 


A/R . (2-1) 

where R is the nose radius. This equation shows that the shock 
standoff distance decreases with Mach number. In this formula it is 
assumed that the detached shock wave is a hyperbola that is 
asymptotic to the free stream Mach angle. The equation for the 
coordinates of the shock around a hemisphere-cone [6] is 


. R.A-R 


(2-2) 




Figure 2-1. Illustration of the flow field in front of the missile's nose. 

B. AERODYNAMIC HEATING AND HEAT TRANSFER 
1. The Heat Transfer Process 

A vehicle operating at supersonic conditions is heated by heat 
generated within the boundary layer that surrounds it. Heat is 
produced by the conversion of kinetic energy through flow 
deceleration in the boundary layer, as shown in Fig.2-2. Heat is 
transferred to the body surface by conduction which is represented 


by Fourier's law of heat conduction given by [7] 





iluated 


ent in the boundary layer is eva! 



Figure 2-2. Thermal boundary layer temperature distribution 
near insulated wall /57- 

diffusion by conduction is a complex three-dimensioi 
problem. The solution of this problem requires the formulation of an 
energy balance which considers the thermal energy stored at each 




well the therm 


energy conducted. 





Figure 2-4. Distribution of temperature near 
Figure 2-3. Effect of heating and cooling on insulated wall [8], 

temperature distribution [8], 

In practice heat transfer is computed using finite differences and 
finite element methods for heat conduction. At elevated temperatures 
a surface will lose heat to its surroundings by radiation according 
to the Stefan-Boltzmann law, given by 

‘iraa ■ (2-5) 

where cs=5.67x10'® Wm'^ K'" and e is the surface emissivity having 
values between 0 and 1 depending on the material, Tbq is the 
background temperature at which the background radiates back to the 
surface. There is also the possibility of radiating from the skin of 
the vehicle to its own internal structure. While the form of the 
Stefan-Boltzmann law is the same, the temperatures of the radiating 
surfaces and non-gray absorption and emission must be considered. 
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Figure 2-5. Heat transfer modes. 


2. Recovery Temperature 

It is important to note that in gases the adiabatic wall 
temperature . is always less than the free-stream stagnation 
temperature, T^. This condition is in contrast to what one may 
expect from the well-known temperature-Mach number relation 


l.[ (y.l)/2]M^ 


(2-6) 


i.[(y-i)/2iM^ 


viscous effects within the boundary layer decrease the velocity 
from the freestream value at the outer edge of the boundary layer to 
zero at the surface. Pure compression heating would give a 
stagnation temperature of 



(2-7) 


the 


rface (that is, in Eq. 2-6 let Tj 


and T2 = T„. Then 











M]^ = 0 and M2 = MJi . is the value of temperature where the 

velocity is zero, i.e., at the stagnation point. On the other hand, 
in continuum flow the velocity at the wall is zero, but as already 
stated, the adiabatic wall temperature is always less than the free- 
stream stagnation temperature 



Figure 2-6. Stagnation temperature versus Mach number at 
several altitudes[9J. 

In practice, heat is conducted outwards through the boundary 
layer. Viscous dissipation produces additional heat and the actual 
temperature reached at the surface (assuming no heat transfer into 
the surface) is called the recovery teraperature[8]. Therefore, for 
practical computations, it is convenient to define a quantity known 
as the recovery factor, r, which is a measure of the fraction of the 
local free-stream dynamic-temperature rise recovered at the wall: 


10 
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It is obvious that if the entire plate is adiabatic, energy 
considerations dictate that the mean stagnation temperature with 
respect to mass flow must be equal to the free-stream stagnation 
temperature Since the wall temperature is less than T^, it 

follows that the distribution of stagnation temperature within the 
boundary layer must be of the form shown in Fig. 2-4, with some 
portion of the boundary layer having a stagnation temperature 
greater than T^. 

The value of the recovery factor depends upon Prandtl number Pr 
which is a ratio measuring the effects of conduction and dissipation 
mentioned above. The Prandtl number is defined by 

Pr . —^ = 0. 71 for air at N.P.T. (2-11) 

where ^ is the coefficient of viscosity and k is the thermal 
conductivity. These coefficients are strongly temperature dependent 
and graphs showing this dependence are drawn in Fig.2-7. It can be 
shown that the recovery factor is approximated in terms of the 
Prandtl number: 

r . Pr " (2-12) 

where n=l/2 for laminar flow and n=l/3 for turbulent flow. While 
this is a useful approximation, the recovery is also a function of 
the pressure gradient and local Mach number. 

The adiabatic wall temperature is then given by an equation such 
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(2-13) 



where is the velocity at the edge of the boundary layer and is 
the specific heat. 

In missile flight the dome is subjected to aerodynamic heating 
resulting from the convective heat-transfer process described by the 
equation 


h . 


(2-14) 


where T„ is the wall temperature and the recovery (adiabatic) wall 
temperature. The adiabatic wall temperature is that temperature at 
which the heat-transfer to the wall is zero, i.e., where the 
temperature gradient in the boundary wall is zero, as shown in 
Fig.2-3. The heat-transfer coefficient is the ratio of the heating 
rate to the thermal driving potential. 

The heat transfer coefficient h in Eq.(2-14) varies considerably 
with flight conditions. It is often more convenient to work with a 
non-dimensional heat transfer coefficient called the Stanton number 
defined by: 


St ■ 




(2-15) 




Stanton number and skin friction coefficient are linked by 


the Reynolds analogy 


(2-16) 


where Cf is the skin-friction coefficient and S is the analogy 
factor. 

3 . Heat Transfer to a Hemisphere 

The aerodynamic flow over a hemispherical nose is laminar at 
the stagnation point. As the air accelerates in passing from the 
stagnation point to the base of the hemisphere, it can go through 
transition, and the boundary layer generally becomes turbulent. 
This, happens at Reynolds numbers around 4x10^ [2] . 

For the laminar case, Sibulkin has shown that locally the 
Stanton number on a hemisphere stagnation point can be written 



where D is the diameter of the hemisphere, 0 is the velocity 
gradient and u is the viscosity. 

All of the fluid properties are evaluated at the local 
conditions just outside the boimdary layer. In the case of 
stagnation this corresponds to stagnation temperature and Pitot 
pressure. This reduces Eq.{2-17) to[2] 

14 



St ■ 


9.37 


(2-18) 


/uTd [ mks 


where (c^ij/kji ^ ^ is the stagnation Prandtl number (approximately 
0.72 for normal flight conditions), p„uD/iJ^ is the free-stream 
Reynolds number based on hemisphere diameter. The subscript o refers 
to total or stagnation values. 

From Ref [1], 



(2-19) 


The subscript 0 denotes the total or stagnation value after the 
air passes through the shock wave. If it is assumed that jj oc 
the perfect gas eguations are used then Eg.(2-18) reduces to the 
general laminar Stanton number equation. This gives the heat- 
transfer rate at the stagnation point where the recovery factor, r. 


h^^/2R - 1. 7Cpy^p.p.a_M. (1.0.2M^) ^ (2-20) 

equals 1. The laminar heat-transfer rate decreases away from the 
stagnation point. At the stagnation point another semi-empirical 
formula can be used to calculate the heat transfer[4] 

From this formula C. Klein showed that the most severe thermal 
environment in terms of temperatures does not correlate with the 
highest thermal-flux situation[4]. 

When the boundary layer is turbulent, the Stanton number has a 
maximum at about 40° from the stagnation point, then decreases 
15 




further aft. 


As shown from the Eq.(2-7) at supersonic speeds the temperatures 
reached are quite high, which raises the issue of suitability of 
certain materials as IR transmissive materials in a high tenperature 
environment. However, it must be determined first whether the time 
frames are such that the surface temperatures can approach 
stagnation-type limits. For that purpose, one may consider the 
dimensionless Biot number 


Si 


hi 


which characterizes the transfer of heat across the boundary layer 
in relation to the dome's ability to transport heat. Here L 
represents the dome thickness. If the wall is thermodynamically 
thin, i.e., Bi«l and the radiation and conduction losses are 
negligible, then the aerodynamic heat-transfer rate can be equated 
to the rate of enthalpy increase within the wall, or 


h(T^-TJ . P„i>„Cp 


(2-22) 


where p„ represents the window density and Cp„ the heat capacity of 
the window. 

If there is a step change in the aerodynamic conditions, e.g. 
velocity (thus a step change in h and T^) , the solution is 



er window temperatures), the radiation heat tram 
lifleant part. If this is to be considered, the h' 
tion becomes 

dr 

-Tj ■ . ae(T)T, . (2 


) can be solved for a series of constant temperature 
e accurate temperature history. The true equilib: 

can also be found from the same equati< 


( dT/dt=0 ) 




Eq.(2-25) is used to give an indication of the relative rate of 


temperature increase. The equation for a time constant, x, is 


P,c^L 

and gives the time to reach equilibrium 


(2-25) 


T^-T I 


0.368 


(2-26) 


in the absence of radiation. In Fig.2-8 the generalized time 
constant for an arbitrary set of conditions, and a typical material, 
MgFj, 0.07 in. thick, is used to get a specific time constant. For 
1 mm thick diamond domes, the heat transfer coefficient at the 
stagnation point leads to very small Biot numbers, which in turn 
suggests thermal time constants of the order of 1 sec at very high 
altitudes. In these conditions, for flight times of about 10 
seconds or more, diamond missile domes will indeed reach the 
recovery temperature[4]. 


C. AEROSTRUCTURE OF THE INFRARED DOME 

As will be discussed in later sections, the self emission of 
the dome window, due to the high teiiperatures induced by aerokinetic 
heating, is an aero-optical effect that will have a major influence 
on the sensor performance in a high speed flight situation. Dome 
emission is going to depend on the emissivity of the dome material 
among other variables. An accurate prediction of the dome emissivity 
18 



should consider the thickness of the dome. Therefore a discussion of 
the calculation of IR dome thickness is relevant. The dome's ability 
to withstand the thermal shock induced by transient heating on a 
fly-out trajectory is one of the major problems encountered in the 
design of high supersonic IR seeker domes. The subject of thermal 
shock resistance of IR domes has been very well characterized and 
analysed by Dr. Claude Klein. The following sub-sections summarize 
some of his published work on the subject[4][10][25]. 

1. Thermal Shock Resistance of the IR Dome 

One of the most constraining factors in the design of IR 
missile window domes is the one attributed to brittle fracture 
caused by thermal shock[10]; the fracture is induced by tensile 
stresses on the "cold" side corresponding to the inner surface of 
the dome. These stresses are caused by the sudden temperature 
gradients generated by aerodynamic heating, which in conjunction 
with thermal expansion, produce compressions on the outer surface 
and tensions on the inside. 

Starting with the expression for the maximum stress experienced 
by a clamped plate, or a complete sphere, 




CLE 

2 <l-\>) 


AT , 


(2-27) 


if there is a linear temperature variation across the thickness and 
if both are free to expand: 


AT ■ \T^^-TJ 


(2-28) 


where AT refers to the temperature difference between the outer and 
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where AT refers to the temperature difference between the outer and 
inner surfaces. In a transient situation, AT is a strong function of 
time and depends critically upon the heat-transport process as 
characterized by the Biot number, Bi. If the condition Bi<l applies, 
AT peaks when the backface temperature starts to rise, which occurs 
long before the front reaches the stagnation temperature, and, thus. 


AT - (T -T.J (Bi) 


(2-29) 


In a first approximation, therefore, a thermally-thin dome will 
be subjected to a peak stress that amounts to[10] 


- 2 ,1 . (2-30) 

where (Bi)^t refers to the Biot number at the stagnation point, T^^ 
is the relevant stagnation temperature, and must be set equal to 
the wall temperature at the onset of the shock. Assessing the 
thermal shock resistance (TSR) then boils down to comparing la^^lp 
with the allowable flexural stress, o>/4, which is best done by 
defining TSR as follows: 


TSR 


10 log 



(2-31) 


if the results are to be quantified. The stress ratio can be 
expressed as a product of three terms. 
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Ap • P^t'P- ■ P-^ (Ps/P.> -^J > 


(2-33) 


where p^c/p^ refers to the pressure ratio across the bow shock and 
along the stagnation streamline: 



The external pressures give rise to a complex stress pattern 
that reflects the nature of the edge attachment, i.e., the 
stiffness of the dome interface with the support ring. Two limiting 
cases will be considered[10] ; (a) the simply supported shell that is 
free to bend when subjected to the action of a uniform normal 
pressure, and (b) the shell with built-in (or clamped) edges that is 
restrained from moving in the "radial" direction. For simply 
supported shells, the maximum stress induced by bending, is 

tensile and obeys the equation 


provided the dome aspect ratio, L/R, satisfies the condition 


sin^(6) ^ sin^ (9) 

12 R 1.2 ' 


(2-36) 


where 6 designates the hemisphere truncation angle (see Fig.2-9). 
For edge-constrained shells, the predominant stress is compressive 
and localized at the base of the dome, near the interface. It 
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relates 


the aspect ratio of the dome 


very simple relation 


sin^(Q) ^ sin^(Q) 

12 ^ 3 


(2-37) 


if the condition 

Ap 


(2-38) 


is met. Evaluating Eq.(2-35) for several angles 8 shows that, on a 
relative scale, the bending stress depends solely on the wall 
thickness (the stress increases as the thickness decreases) and can 
be approximated by means of the relation 



(2-39) 


as long as the conditions (2-37) are satisfied. It is also seen that 
this stress represents a worst case in the sense that constraining 
the edge displacement would decrease the stresses, especially for 
low aspect ratio domes. Because of the thermal expansion mismatch at 
the interface of the infrared dome and its support structure, a 
"floating" attachment scheme that decouples the dome from the 
missile body would be desirable, and such optimized attachment 
should result in pressure-induced tensile stresses as given by 
Eq.(2-39). The maximum tensile stresses failure criterion, i.e.. 


where o> refers 1 


I the nominal flexure strength of the dome material. 
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and SF is an appropriate safety factor, then yields the minimum 
thickness required to prevent fracture. Since currently available IR 
transmitting dome materials are brittle ceramics, they may exhibit 
a wide range of fracture probabilities under apparently identical 
loadings, which suggests safety factors of at least four (4) if 
is the characteristic flexural strength derived from Weibull 
statistical analysis. It is concluded that the dome thickness 
required for withstanding a pressure load Ap with good probability 
of survival can be obtained from the following equation: 



D. RADIATION PROPERTIES OF SURFACES 

In analysing the interchange of radiant energy among surfaces, 
it is necessary to know their emission, reflection, and absorption 
characteristics. Besides the obvious need for transparency, a 
critical requirement for an IR window in a hot environment is that 
it emits little radiation that would obscure the scene being viewed. 
In the next sections the emission properties of windows will be 
explored. 

1. Radiometry 

In describing the directional transport of energy, it is 
convenient to deal with the radiant energy (in joules) contained 
within a small solid angle dcj. This naturally leads to the concept 
of the radiant intensity. The radiant flux, 0 (in watts) emitted 
per unit solid angle(w) by a point source in a given direction is 
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called radiant intensity, I, and is given bY[19] 


d3> 

do) 


[W/sr] . 


(2-42) 


Another useful radiometric term is the radiant flux density at a 
surface, measured in watts per square meter. It may be either 
emitted from a surface , in which case it is called radiant 
exitance, M, or incident onto a surface, in which case is called 
irradiance, E. 

The inverse-square law of radiation from a point source is 
derived by calculating the irradiance of a point source on a 
spherical surface surrounding the point, of solid angle An sr and 
surface area 4nr^. Hence 


^ 4nJ 
A 4nr^ 


(2-43) 


The radiance, L, represents the radiant energy per unit time 
and unit surface area that passes through a solid angle da inclined 
at an angle 9 with respect to the surface normal. This situation is 
represented in Fig. 2-10. The radiance is given by 


dl 

(dA cosQ) 


[N/m^-sr] 


(2-44) 


In the case of spectral quantities (wavelength-dependent 
distribution), a subscript A is affixed to the radiometric 
quantities, whereas for total radiation these subscripts are 
omitted. In general, the direction of the intensity is specified by 
angles 6 and 0 of a spherical surface as shown in Fig. 2-10. In the 
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case of a Lambertian surface, the intensity for a certain radius of 
the surface varies in accordance with the following relation, 

1(d) . 1(0) cos e . (2-45) 

In this case the radiance is independent of the viewing angle, and 
the surface is said to be perfectly diffuse. The intensity 
distributions are related to the corresponding hemispherical 
radiation fluxes by integration over the hemisphere surrounding the 
surface. As an example, for a flat surface, the radiant power 0 
leaving the surface is connected with the radiant intensity at the 
surface as follows: 

<!>. / 1 ( 0 ) cos e dm [W] 

in which the subscript 'hem' denotes the hemisphere. Similar 
relationships apply for the incident and reflected radiant fluxes. 


n 



Figure 2-10. Angle Notation. 
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2. Blackbody Radiation 

The ideal blackbody is a perfectly efficient and diffuse 
emitter of radiant energy. No body at the same temperature can emit 
more radiation at any wavelength or into any direction than can a 
blackbody at the same temperature. It follows that a blackbody is 
also a perfect absorber: all radiation incident on a blackbody will 
be absorbed regardless of wavelength or angle of incidence. The 
spectral radiant exitance of a blackbody is given by the Planck 
equation, 


- 


2nhc^ 1 

AS j. 


(2-47) 


where the physical constants h, c, and k represent the Planck 
constant, the speed of light, and the Boltzmann constant 
respectively. After substituting the values of these constants in 
the Planck equation, the result is 


3. 745 X 10^ 


(2-48) 


where A is in micrometers and T is in kelvins. The spectral radiant 
exitance for three background temperatures is plotted in Fig. 2-11. 

The peak exitance, shifts toward shorter wavelengths with 

increasing temperature according to the Wien displacement law, 
expressed as 

^ ^ [vun-.ff] ■ 

Integration of the spectral radiant exitance of Eq.(2-47) over all 
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wavelengths gives the total radiant 


tance given by 


(2-50) 


known as the Stefan-Boltzmann law, with a as the Stefan-Boltzmann 
constant, equal to 5.67xlO'®W/m^-K''. 

The emissivity z expresses quantitatively the fact that real 
surfaces always emit less radiation then a blackbody. The blackbody 
emissivity is expressed as the ratio of the total radiant exitance 
of the body to the total radiant exitance of a blackbody at the same 
temperature. 


e (T) . 


M 


(2-51) 
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Figure 2-11. Blackbody radiation spectral distribution at three background 
temperatures [20]. 

Directional Emissivity 

ingular emissivity is defined as the ratio of the 
Lcular direction to the radiance of a blackbodi 
:tion. In general, for arbitrary surface condit; 
.on of the emitted radiance L depends on two ang 


Fig. 2-10. On 


nd, the blackbody 





Lgg is uniform. Consequently, one can write 

e . (2-52) 

Eg. (2-52) states that, for a diffuse surface, the directional 
emissivity is independent of angle and is equal to the hemispherical 
emissivity. 

The directional emissivity characteristics of real surfaces are 
revealed by measurements performed on a variety of materials by E. 
Schmidt and E. Eckert[11]. One of their measured distributions is 
reproduced in Fig. 2-11. The figures are in the form of polar 
diagrams in which the concentric circles are contours of constant 
emissivity s and the radial lines are contours of constant angle 6. 
The surfaces employed in the experiments were presumably isotropic, 
and no dependence on the angle $ is indicated. The temperatures of 
the various test surfaces were typically a few hundred degrees 
Fahrenheit. The directional emissivity results were determined from 
total radiation measurements. 

Figure 2-11 shows that, for nonmetals, the directional 
emissivity remains essentially uniform for inclination angles 
between 0° and 50-60°; then it drops off quite rapidly to zero. In 
the region of near uniformity, the emissivity values are in the 
range from 0.8 to 0.95, which is typical for dielectrics. 
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Figure 2-12. Distribution of the total directional emittance for several dielectrics[11 ]. 

The directional eitiissivity characteristics displayed by 
metallic surfaces are somewhat different from those just discussed. 
In the range of angles between 0° and 30-40°, e remains quite 
uniform. However, at larger angular inclinations e increases quite 
sharply. The percentage increase is more pronounced when the level 
of emissivity is low. At angles approaching 90° the directional 
emissivity should approach zero. 

The directional absorptivity is the fraction of the radiant 
energy incident on the surface from a specific direction that is 
absorbed. In accordance with Kirchoff's law, the directional 
monochromatic absorptivity and the directional monochromatic 
emissivity are equal provided that the incident beam is uniformly 
polarized. That is, 

o(^ (e,(b) • ce,<t); . (2-53) 


3J 







4. Optical Properties of Dielectric Window Materials 
In the infrared, there are two main contributors to the complex 
index of refraction of insulating ceramics. They are electronic 
transitions and lattice vibrations or phonons (a quantized packet of 
vibrational energy)[12]. The electronic transitions contribute to 
the real part of the index of refraction only in the form of a 
frequency-independent parameter with weak temperature dependence. 
Most of the frequency and temperature dependence of the complex 
index of refraction at mid-infrared frequencies comes from phonons. 

Single-phonon absorption processes are due to the excitation of 
a phonon upon the absorption of a photon. These processes are very 
strong, since the photons can couple directly to the fundamental 
lattice vibrations of a material. Multiphonon absorption processes 
are due the excitation of more than one phonon upon absorption of a 
photon. These processes become important in the transparent regions 
of a material where one-phonon absorption processes do not dominate. 
Difference bands of multiphonon absorption occur in the millimeter 
spectral region, and sum bands occur in the infrared region. Sum 
band processes typically dominate infrared emissivity[12]. 

The behavior of the absorption coefficient and refractive index 
of a typical optical material is indicated schematically in Fig. 2- 
12. In general, materials absorb light strongly in the infrared and 
ultraviolet regions, near wavenumbers of 10^ and 10^ cm‘^, 
respectively, in the upper part of Fig. 2-13 [12] . Infrared 
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Figure 2-13. Absorption coefficient and refractive index of a typical optical window material. Adapted 
from Thomas[9]. 

radiation stimulates vibrations of the atoms in a solid at the 
characteristic Reststrahlen frequencies. Ultraviolet radiation 
promotes electrons from filled bands to higher energy unfilled 
bands. Between the two strong absorption regions is a window in 
which the material is transparent. The infrared edge of the window 
is called the multiphonon region because the transitions responsible 
for the weak absorption arise from multiple excitations of one or 
more vibrational modes(phonons). The ultraviolet edge of the window 
is called the Urbach tail. 

The optical properties of non-magnetic dielectric window 
materials can be specified completely by the complex index of 
refraction, n(v,T) or complex permittivity, 6’(v,T) where v is the 
wavenumber (reciprocal wavelength) and T is temperature[12]. The 
relationship between n and e is 
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(\),T) - ik(V,T) 


(2-54) 


n(V,T) 


where 

e(v; - (n^ - k^) - 2ink ■ - ie^ , (2-55) 

n is the index of refraction, and k is the index of extinction. 

The magnitude of the reflectance of the medium (at normal 
incidence) is 


The absorption coefficient is defined by 


I (v,r; - 4nk(v,T) . 


The radiance from a material, L„ (v,T), is expressed as 
L^(V,T) . £ (e,V,T)L^^(V,T) 


where s(9,v,T) is the spectral emissivity of the material, Lgg(v,T) 
is the radiance of a blackbody, T is the temperature, v is the 
spectral,frequency and 6 is the angle of incidence (i.e. the angle 
between the surface normal and the incident ray). In general, the 
emissivity is composed of two components: a specular and a diffuse 
component. For our purposes, the diffuse component is negligible 
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since the materials have smooth flat surfaces and negligible bulk 
scattering. 

Based on the Stefan-Boltzmann law and Kirchof f' s law of 
radiation, the specular component of emissivity can be expressed 
as[12] 




T) 


(1 - R(Q,V,T)) (1 . 
1 . 


(2-59) 


where R is the single-surface Fresnel power reflectance for 
unpolarized light, and t is the window thickness. In general £•„ and 
R have angular dependencies consisting of sine and cosine functions. 
Near normal incidence (i.e. near 6=0° ), the angular dependencies of 
and R can be neglected, since the sine and cosine are 
approximately constant. 

Usually, a measurement of |i?l is not sufficient to characterise 
£. Therefore, it is convenient to describe the form of the complex 
permittivity. The classical pole-fit model represents the 
fundamental lattice vibrations or single-phonon transitions in a 
material. This feature is important in representing the index of 
refraction at all optical frequencies, but the absorption 
coefficient obtained from the model is valid only in the vicinity of 
the one-phonon frequencies. The absorption coefficient needs to be 
known in the region of transparency as well, requiring a description 
of multiphonon absorption bands in which more than one-phonon 
excitation occurs upon the absorption of a single photon. 




E. DOME EMISSION AND EMISSIVITY AT HIGH TEMPERATURES 

As the IR dome heats up, it will emit graybody radiation over 
the entire electromagnetic spectrum, which may adversely impact the 
operation of the seeker system in a high-speed flight situation 
[2] [4] [13]. It is recognized that saturation of the detectors by 
radiation from aerodynamically heated domes is not a major issue but 
the added background-induced noise will lower the signal-to-noise 
ratio for detection, decreasing the effective performance of 
the system in terms of detection range as well as tracking ability. 
An accurate prediction of the dome emission requires a good 
knowledge of the dome emissivity properties at the high operating 
temperature. The emissivity can be obtained by models or by 
experimental data. 

The effect of temperature on the emissivity properties of 
several dome materials is shown in Fig. 2-14. The absorption edge 
of an infrared window material shifts to a shorter wavelength as the 
temperature increases (Fig. 2-14). This occurs as the weak 
absorptions present in this region begin to dominate. 

When temperature of a material increases(see Fig. 2-15), so does 
its emissivity. Fig. 2-16 shows the effect predicted for 
sapphire [9] . As shown in Fig. 2-15 the sapphire window at 1000°C 
(1273 K) emits 7.4 times as much radiant energy as the same window 
at 500°C (773 K) if the emissivity is independent of T[9] . Near the 
absorption edge the emissivity increases with temperature, so the 
actual emission is even more than 7.4 times the lower temperature 


36 




Wavelength (;im) 


Figure 2-14. Emissivily of : 
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Figure 2-1 7 . Emissivity of 0.20 cm thick sapphire calculated from 
measured absorption[9]. 



Figure 2-18. High-temperature spectral emissivity of a hemex grade sapphire o-ray sample (12 .7 mm 
diameter, 1.0 mm thick) over the temperature range of950 K to 1930K. [12], 
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111. GOVERNING EQUATIONS OF FLUID MECHANICS 


In order to compute compressible viscous flow around a body, 
the continuity, momentum and energy equations must be solved 
simultaneously. Here the vector and the conservative-law form of 
the compressible Reynolds-averaged Navier-Stokes equation is 
presented. 

A. CONTINUITY EQUATION 

The continuity equation expresses the Conservation of Mass law. 
This law applied to a fluid passing through an infinitesimal control 
volume fixed in space yields the following equation of continuity 

^ . V ■( pv ). 0 . (3-1) 

Here p is the fluid density and V is the fluid velocity vector. 
Eq. (3-1) states that the net mass flux through a control volume 
surface must be equal to the time rate of change inside the control 
volume. For a three-dimensional Cartesian system, where u, v, w 
represent the x ,y, z components of the velocity vector, Eq. (3-1) 
becomes 

3p d d d 

Jl ■ ^ ° • (3-2) 


4J 




B. MOMENTUM EQUATION 

The momentum equation expresses Newton's second law as applied 
to a fluid element passing through an infinitesimal, fixed control 
volume in space. The momentum equation is: 

C pv ; . V • pw . pf . V • . (3.3) 

The first term in Eq.(3-2) represents the time rate of change 
of momentum per unit volume in the control volume. The second term 
is a vector representing the momentum convection, per unit volume, 
through the control surface. The first term on the right-hand side 
is the body force per unit volume and Ilij is the stress tensor given 
in compact tensor notation by 

( au. au 1 au„ 

where 5^^ is the Kronecker delta symbol (5^^ = 1 if i = j and = 0 
i.f i ^ j ) ; , U2 , Uj represent the three components of the 
velocity vector V; u is the coefficient of viscosity (dynamic 
viscosity) and X is the second coefficient of viscosity. 

By substituting Eq.(3-4) into Eq.(3-3) and writing the Cartesian 
component velocities in the directions x, y, z as u, v and w, for 
clarity, the Navier-Stokes equation is obtained. In component form 
the Navier-Stokes equations are; 
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C. ENERGY EQUATION 

The energy equation is obtained by applying the First Law of 
Thermodynamics to a fluid passing through an infinitesimal, fixed 
control volume 


1£ 


. pf • V . V • ('ll ■ v; 


where e,. is the total energy per unit volume. The first term in left 
hand-side of Eg.(3-8) represents the rate increase of total energy 
per unit volume in the control volume while the second term 
represents the rate of total energy flux per unit volume through the 
control surface. The first term on the right hand-side is the rate 
of heat produced by external agents where Q is the heat added to 
the system. Vg is the rate of heat conduction through the surface. 
Fourier's law for heat transfer by conduction states that the heat 
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transfer q 


be expressed as 


g = -kVT (3-9) 

where k is the coefficient of thermal conductivity and T is the 
temperature, f represents the body forces per unit mass, pf-V is 
the work done on the control volume by body surfaces, and V ~V) 
is the work done on the control volume by surface forces. 

The conservation of energy for a three-dimensional Cartesian 
coordinate system becomes 

- p (f^u + fyV + f^w) + (e^u + pu ~ q 

* pv - - VV^ + + gj 

dy 

+ -5- + V pw - ur _ - vr -WT_, + qj = 0 

oz ^ 

(3-10) 


where 


-i-f) 

-i) 


(3-11) 



D. EQUATIONS IN NONDIMENSIONAL VECTOR FORM 

In order to facilitate the numerical implementation the 
governing fluid dynamic equations are written for a body fitted 
coordinate system. The equations are transformed to a generalized 
body fitted coordinate system using the following transformations 

5- \ {x,y,z), r| . r| rx,y,z; , l,{x,y,z) (3-12) 



PHYSICAL DOMAIN COMPUTATIONAL DOMAIN 



Figure 3-1. Transformations and body coordinate system[15]. 


The strong conservation law is retained for a curvilinear 
coordinate system p , ^ and the governing equations become[15] 
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d^{E . e\) . d^{F . F^) . (G . G^) - 0 


(3-13) 
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(e^.p) a ■ 


(e^.p)V. n^J 


(e^.p)H - 4^ 


U- 

V. n^.n^u.n^v.n^., O.is) 


where U, V, and W are contravariant velocities written without 
metric normalization. 

The viscous terms are given by 








' <, • J-' 
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- X (u^.v^.w J ■2}JU^ 

■ fj (u^.V^) 

T„ • . p ■ V,J 

Tyy - X (u^ ■ . wj . 2]2V^ 

■ '^zy ■ ^ ^^z ■ 

- X Cu^. v^. w^; . 2pw^ 

3^ . YK Pr'^ d^e^ . 

3^ - yk Pr'^ a^e^. . 

3^ ■ YK Pr'^ d^e^ . 
e^ . e^p-^ - ^ Cu^ . 


(3-17) 


Here Cartesian derivatives are expanded in terms of ^ , p , ^ 
coordinates by the chain rule relations as 


(3-18) 


The Cartesian velocity components u, v, w are non- 
dimensionalized with respect to a« (the free-stream speed of sound), 
density, p, is referenced to p«. and total energy, e, to 
p^aj. Pressure is given by 

P ■ (Y-l) . I p (u^v^.w^)] (3.19) 

and throughout y is the ratio of specific heats. Also, k is the 
coefficient of thermal conductivity, p is the dynamic viscosity, 
while A from the Stokes' hypothesis is -(2/3)u. 
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The metric terms and the Jacobian of the transformation are 
obtained from chain rule expansion of , etc., and solved for 
ix ’ iy ' etc. These terms are given by. 


CynZ^-y^zii; 

^x 



by 





^x ■ 

5, 







■ -^x^x-y.^y-^T?z 





E. THIN LAYER APPROXIMATION 

In order to accurately compute the normal gradients in the 
boundary layer, it is necessary to make the normal grid spacing very 
fine close to the body surface. In the streamwise direction the flow 
gradients are not large, when compared to those in the normal 
direction, and the resolution of the grid can be lower. This results 
in a mesh where the cells have a high aspect ratio close to the body 
and a smaller aspect ratio when approaching the outer boundary of 
the computational domain. This procedure is going to save 
considerable computer storage and computation time without 
sacrificing the accuracy of the computation. 

In a thin-layer approximation all viscous derivatives in the E, 
and r direction (along the body) are neglected, while terms in 
normal direction ^ are retained, (see Fig.3-1). Eq. (3-1) thus 
simplifies to 
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(3-22) 


B^q-. B^e' . B^F. 




Re-‘B^ S~ 


0 

l^(C- Vy Cz) (u/3) 

U (Cy Vr) (11/3) 
q (iy ll- O «c (P/3) 

^l> [0.511 (u/v^.w^) ^.KPr ' (Y-l)'^ (a^) ,_} 

. (^,/3) ((n/3) (^^u^.i:^v^.[,w^) \ 


Unlike the boundary-layer model the thin layer accounts for 
pressure variations in the normal direction and does not require 
matching an inviscid layer with a viscous layer. 

The eddy viscosities were computed using the algebraic Baldwin- 
Lomax eddy viscosity model. 
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IV. FLOW FIELD COMPUTATION 


The flow field was computated in order to obtain the temperature 
variation on the surface of the dome, corresponding to the IR 
window. The aerodynamic heating of the infrared window is 
responsible for a degradation of the seeker infrared sensor 
performance at supersonic speeds. 

A. METHODOLOGY 

For the space discretization the Steger-Warming flux vector 
splitting is used along the streamwise direction. Along the 
circumferential and normal directions central differences with added 
numerical dissipation are used. The integration is performed with a 
two-factor partially flux split implicit algorithm. A combination of 
second and fourth order differences have been used for smoothing 
along the central-differencing directions. Fourth-order dissipation 
terms are added explicitly to the equations to damp high frequency 
growth and thus serve to control nonlinear instability. The 
resulting implicit code is second-order accurate in space. 

A body-fitted grid, clustered in the radial direction was 
employed to resolve the viscous layer adjacent to the body surface. 
This grid was obtained by means of a hyperbolic grid generator. 

The flow field was computed using the OVERFLOW Version 1.6t 
code developed at NASA Ames Research Center running in a STARDENT 
workstation installed at the Computational Fluid Dynamics Laboratory 
(Naval Postgraduate School). After a solution was obtained, the 
original grid points were redistributed (adapted) based on the 
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initial flowfield solution, and a new solution was computed using 
the adapted grid. This procedure increases the accuracy of the flow 
prediction by decreasing the effects of flowfield discretization. 
This procedure is particularly important, in supersonic flow 
computations, where the existence of strong gradients and shocks 
require careful distribution of grid points to minimize errors and 
produce more accurate flow field predictions. 

B. BODY GEOMETRY AND GRID GENERATION 

The flowfield was computed around three distinct body 
geometries, namely two nose shapes (cases I & II) and a test 
geometry (case III) . The case III geometry had the shape of a 
hemisphere-cylinder. The computed shock stand-off distances and 
surface pressure distribution of this latter configuration was 
compared with experimental data to give confidence in the 
computations. 

Several C-0 type grids were generated, corresponding to the 
hemispherical-ogive (case I) and to the ellipsoidal-ogive (case II) 
forebody. Longitudinal cross sections of the computational grids for 
the case I and case II are shown in Fig. 4-3 and Fig. 4-4 
respectively. The case I body consists of a truncated hemisphere 
forebody(1.9024 cm long) attached to a 3:1 tangent ogive (12.728 cm 
long) attached to a 4 caliber tangent cylindrical afterbody. The 
case II body consists of a truncated ellipsoidal forebody(2.6943 cm 
long) attached to a 3:1 tangent ogive (12.72 8 cm long) 
attached to a 4 caliber tangent cylindrical afterbody. The 
hemispherical and ellipsoidal portions of the nose correspond to the 



IR window locations respectively. The origin of the physical 
coordinates is at the nose tip for both cases. The three- 
dimensional body fitted grids were generated using the HYPGEN code. 
After a flow solution was found an adaptative grid code[17], SAGE, 
based on the algebraic, solution-adaptive scheme of Nakahashi and 
Deiwert, was used, in order to enhance the resolution for better 
capture of the shock front. A third inviscid grid was designed to 
compute the flow around a hemisphere-cylinder(case III).- 
corresponding to the test case. The geometrical shapes of the bodies 
considered are as follows, 

1. Case I. Hemisphere-Ogive-Cylinder (Fig. 4-1) 

This nose configuration corresponding to the hemispherical dome 
is defined using the following relation 

r . \/3.3166^- (z-3.3166)^ , Oszil.9024 

r. - (z.3.3697)- ^ 1.9024szsl4.6303 (4-1) 

3 54 

r . 6.0000 , Zil4.6303 

The origin of the coordinate system is located at the tip of 
the nose. Fig. 4-3 shows the computational grid used for the 
calculations. The resolution of the grid is increased around the 
shock wave where the flow gradients are large. 
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Ellipsoidal Window anached to a Parabolic Nose Shape 



z, axial distance 

Figure 4-1. Hemisphere-Ogive nose configuration. 


Ellipsoidal Window attached to a Parabolic Nose Shape 



z, axial distance 

Figure 4-2. Ellipsoid-Ogive nose configuration. 
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II. Ellipsoid-Ogive-Cylinder 


2. Case 

This nose configuration corresponding to the ellipsoidal dome 
is defined using the following relations: 
elipsoid excentricity = 0.9 


. - (z.2.577)-- 
. 6.0000 , 


- , 2.6943<:zsl5.4223 
Z215.4223 . 


Figure 4-4 shows the 77x93x31 computational grid used for the 
calculations. 

3. Case III. Hemisphere-Cylinder (Fig. 4-2) 

This body configuration corresponds to the test case. Its shape 
is defined by the following relations: 


r . ^36.00^- (z-6.00) ^ 
r . 6.00 , 


C. GRID ADAPTATION 

Grid adaptation procedures in general either redistribute 
points or refine distribution (by adding points) to reduce solution 
errors by minimizing grid discretization errors. Gnoffo[17] first 
introduced a method for the redistribution of grid points based on 
local flow gradients. This method is analogous to finding the 
equilibrium position of a system of springs that connect each node 
with adjacent nodes and whose spring forces (or tension constant) 







Figure 4-3. 77x48x8! viscous computalional adapted grid for the Hemisphere-Ogive-Cylinder. 
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Figure 4-4. Three dimensional Hemisphere-Ogive-Cylinder computational mesh. 
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is proportional to the local error or weight function. The 
multidimensional problem of grid adaptation is split into a series 
of one-dimensional (1-D) problems along the computation coordinate 
lines. The reduced 1-D problem then requires a tridiagonal solver 
to find the location of grid points along a given coordinate line. 
Multidirectional adaptation is achieved by the sequential 
application of the method in each coordinate direction (Fig. 4-5) . 



Figure 4-5. Splitting of2-D adaptation into two 1-directional adaptations[17], 


In the adaptive grid procedure outlined by Nakahashi and Deiwert 
(1985) the tension force directs the redistribution of points to the 
strong gradient regions. To maintain smoothness and a measure of 
orthogonality of grid lines, torsional forces are introduced that 
relate information between the family of lines adjacent to one 
another[16]. The smoothness and orthogonality constraints are 
direction-dependent, since they relate only the coordinate lines 
that are being adapted to the neighboring lines that have already 
58 


I 







been adapted. This implies that the solutions are nonunique and 
depend on the order and direction of adaptation. Nonuniqueness of 
the adapted grid is acceptable since it makes possible an overall 
and local error reduction through grid redistribution. 

D. RESULTS AND DISCUSSION 

Table 4-1 shows the computed values of absolute temperature, 
normalized temperature and coefficient of pressure along the 
stagnation streamline obtained by the solution of the thin-layer 
Navier-Stokes equations about the hemisphere-ogive-cylinder nose 
(Case I) at zero incidence at Mach 4. A 77x7x35 grid point inviscid 
grid was used. The solution was considered converged after 1500 
steps. Figs.4-6 and Fig.4-7 show the plot of these flow quantities 
versus distance from the stagnation point. The bow shock stand-off 
distance taken from the computed flow field is approximately 0.36cm, 
corresponding to a A/R=.ll. From the plot of Fig. 4-7 we can see 
that the temperature jump across the shock is 4.15 compared with the 
value of 4.04 obtained from normal shock tables. 

Table 4-2 shows the computed values of absolute temperature, 
normalized temperature and coefficient of pressure on the surface of 
the nose along a streamline for the hemisphere-ogive at Mach 4 for 
the same conditions of Table 4-1. Fig. 4-8 shows the computed 
temperature profile of the dome surface temperature. The computed 
normalized stagnation temperature is 4.6 comparing with 4.2 obtained 
using Eq.(2-7). Fig. 4-9 shows the computed dome pressure 
distribution. Fig. 4-10 shows the OVERFLOW input file used to 
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generate the values of Tables 4-1 and 4-2. Fig. 4-11 shows the 
three-dimensional flow-field around the missile afterbody. The IR 
dome corresponds to the white portion of the nose tip. Fig. 4-12 
shows the density field in front of the IR dome. 

The flow field has cylindrical symmetry at zero angle-of- 
attack. With missile incidence the flow is no longer axisymmetric. 
Figs. 4-13 and 4-14 show the Mach number distribution and the dome 
temperature field in the dome for an angle of attack of 5°. Fig. 4-15 
gives the density field. It was obtained using the viscous adapted 
grid of Fig.4-5. 

In the computation of the flow it was assumed that the dome had 
reached thermal equilibrium and the wall was adiabatic. In order to 
obtain the equilibrium temperature the dome should be divided into 
elements and a procedure like the one outlined in Eq.(2-24) should 
be carried out to consider radiation and conduction losses. The heat 
transfer coefficient at stagnation is given by Eq.(2-20). For other 
elements that make an off-axis angle (angle between the longitudinal 
axis of the missile and the line between the dome element and the 
dome center of curvature) , d, the value of should be multiplied 
by cos (9)^^^. In this study the dome material was assumed to be 
sapphire or CVD(Chemical Vapour Deposition) Diamond. These materials 
make it possible to fabricate IR domes with a high thermal shock 
resistance(see Appendix A) suitable for the MWIR(Medium Wave 
Infrared) and LWIR(Long Wave Infrared) spectral bands respectively. 
In Appendix A it is shown that for a diamond window at Mach 4 a dome 
thickness of 1mm is enough to withstand the flight conditions. 
Therefore we can expect that conduction losses are not going to be 
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very high since the area of the attachment is small. Also, 
emissivity of these materials is small we can expect small radiation 
losses. In this case, considering the wall adiabatic might be a good 
approximation. 


Table 4-1. Computed streamline flow field properties (Mach-4.0, a=0). 
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265854 
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4.075641 
3.453679 
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001526 
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1.769331 
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Table 4-2. Computed hemispherical dome temperature profile (Mach=4.0, a=0) (adiabatic 
wall). 
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Figure 4-6. Plot of the coefficient of Pressure along the stagnation streamline for the computed flow field 
around the hemisphere-ogive-cylinder 



Figure 4-7. Plot of the normalized temperature along the stagnation streamline for the computed flow jieid 
for the hemisphere-ogive-cylinder. 


















Figure 4-11. Thin-Layer Navier-Stokes Mach Number flow field solution for the Hemisphere-Ogive body 
(Mach=4.0, 77x7x35 computational grid). 
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Figure 4-12. Gradient of index "lens "created by the compressed gas in front of the dome at Mach=4.0. 
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Figure 4-/3. Thin-Layer Navier-Slo/ces Mach Number flow field solution for the Hemisphere-Ogive body i 
5''angle-of-attack (Mach=4.0, 77x48x81 adapted grid). 











Figure 4-14. Hemispherical dome temperature field for a Mach 4 seeker with 5" angle-of-attack obtained 
by solving the Thin-Layer Navier-Stokes equations. 


69 



















Figure 4-15. Thin-Layer Navier-Slokes density flow field solution for theHemisphere-Ogive body at 5° angle- 
of-attack (Mach=4.0, 77x48x81 adapted grid). 
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E. VALIDATION 


In order to validate the computation, a comparison between the 
predicted flow field about the missile nose and experiments was 
tried. After an exhaustive library search no data suitable for 
comparison with the computation was found in the open literature. 
Most experiments were performed at hypersonic speeds. The only 
suitable experimental data to compare with the flow field 
computation were ones for a hemisphere-cylinder at M, = 1.2 and a 
Re = 5.4 X 10® ft. For this comparison a 61x7x41 inviscid grid was 
created. The experimental data was taken from shadowgraphs obtained 
in the AEDC Aerodynamic Wind Tunnel (IT)[18]. 

A plot of the shock standoff distance at both the body axis and 
the wind axis for a hemisphere-cylinder at Mach 1.2 is shown in Fig. 
4-17. The normalized standoff distance, a/R for zero incidence is 
-1.26. The flow was computed using a 61x7x41 inviscid grid. The 
computed normalized standoff distance was -1.0. The poor result is 
considered to be due to lack of resolution of the grid near the 
shock front. Due to restrictions of time another grid was not 
generated to try another comparison. Fig. 4-16 shows the computed 
Mach number distributions using the inviscid Euler solution. 
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Figure 4-16. Hemisphere-Cylinder Euler solution of the Mach Number flow field (61x7x41 grid). 
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-From Shadowgraph 



Figure 4-1 7. Shock standoff distances about Hemisphere-Cylinder at 
incidence from shadowgraphsflSJ. 


73 






V. INFRARED SEEKER DETECTION PERFORMANCE DEGRADATION 

The major effect of the increase in dome temperature on the IR 
seeker sensing performance at high supersonic speeds is considered 
to be thermal emission to the detectors. This increased dome 
radiation contributes to the noise and thus increases the Noise 
Equivalent Irradiance {NEI) to the detector or Noise Equivalent 
Quantum Flux {NEQF) to the detector array, decreasing the effective 
detectivity of the vehicular sensor and thus decreasing the 
detection range of the seeker. 

Initially the basic figure of merit parameters will be 
introduced. Then a discussion of the range equations for non-imaging 
and imaging sensors in a high speed flight environment will follow. 
Several computations of detection ranges for a staring sensor are 
carried out for a range of mean dome temperatures, representing a 
range of missile speeds up to Mach 4. 

A. SENSOR PERFORMANCE PARAMETERS 

It is convenient to introduce the seeker sensor figures of 
merit before starting an analysis of the seeker performance. 

1. Basic Figures of Merit 

A basic figure of merit that applies to all detectors with 
electrical output is responsivity. Responsivity is the ratio of the 
output (usually in amperes or volts) to the radiant input (in 
watts) . The spectral voltage responsivity of a detector at a given 
wavelength A is the measured voltage output, V^, divided by the 
spectral radiant power incident on the detector, Og(A)[29]: 
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(5-1) 


where f is the modulation frequency. 

The noise equivalent power {NEP) of a detector is the required 
power incident on the detector to produce an average signal output 
equal to the rms noise output. Stated another way, the NEP is the 
signal power level that produces a signal-to-noise ratio of 1. The 
voltage signal output is 


so the signal-to-noise ratio is 


(5-2) 


s/w . 


(5-3) 


The NEP is the incident radiant power, <1)^, for a signal-to-noise 
ratio of 1: 


(5-4) 


Solving for NEP gives 


NEP ■ 


(5-5) 


where v,^ is the root-mean-square noise voltage in volts and is 
the voltage responsivity in volts per watt. 
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The detectivity D of a detector 
the noise ecjuivalent power: 


defined as the reciprocal of 


[W'J- (5-6) 

A more useful performance parameter is the normalized detectivity D" 
(pronounced dee-star) which normalizes D to the detector area and 
noise bandwidth, Af, as follows. 


- oJKJf . , 

V ^ NEP 


(5-7) 


given in units of cm-Hz^^^/W. 

The NEP can be divided by the area of the collecting aperture, A^, 
leading to a parameter called noise equivalent irradiance {NBI) [27], 


NET ■ 


NEP 


(5-8) 


2. Noise Equivalent Ten^serature Difference (NETD) 

The most widely used radiometric measure of the ability of 
a military mid-IR system to discriminate small signals in noise is 
the Noise-Equivalent Temperature Difference or Noise-Equivalent 
Differential Temperature. As defined by J. M. Lloyd [22], " The NETD 
is the blackbody target-to-background temperature difference in a 
standard test pattern which produces a peak-signal to rms-noise 
ratio {SNR) of one at the output of a reference electronic filter 
when the system views the test pattern". 

In terms of system parameters the ratio of the differential 
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signal voltage to the detector rms noise voltage in the detector 
be expressed, for a rectangular detector, as[22] 


Av^ 

V 


r 

/ -^ D-(\) T (\) d\ 

J dT ° 


(5-9) 


where a, /3 = detector angular subtenses [radian] 

a, b = detector dimensions [cm] 

D*(A) = specific detectivity as a function of 
wavelength , evaluated at the electrical 
frequency at which the noise voltage spectrum of the 
detector is normalized to unity [cm / W] 

= equivalent noise bandwidth [Hz] of the NETD test 
reference filter 

= effective area of the collecting aperture [cm^] 

Mj^(A) = spectral radiant exitance [W/pm-cm^] 
t^(A) = spectral transmittance of the optics. 

The NETD is defined as the temperature difference for unity SNR, so 
we can set AV^IV^ = 1 and solve for the resulting AT, which is by 
definition the NETD-. 



: B M. 

ap A f - - D-(X) T„(X) dX 

° J BT ° 


(5-10) 


The spectral radiant exitance, , is defined by Planck's 

blackbody radiation law 
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This last equation is more conveniently expressed in • 
[W/cm^ym] for A in [ym] as: 


where = 3.7415 x 10^ [watt ym* / cm^ ] 

Ci = 1.4388 X 10’ [ym K]. 

The noise equivalent temperature difference, NETD, representing 
the total output temporal and pixel(spatial) noise at a given 
background temperature converted into an equivalent temperature 
difference at the background scene, is a very effective figure of 
merit for a staring focal plane array (FPA). 

3 . Noise Equivalent Quantum Flux for a CCD Array 
To develop such an array figure of merit (FOM), an appropriate 
performance measure for a single pixel should be introduced first. 
The figure of merit for a charge couple device (CCD) array should 
take into consideration the filling efficiency of the array, spatial 
resolution limit, modulation transfer function, quantum efficiency 
and integration time[29]. A performance measure for a single pixel 
may be used as a normalization factor for the array modulation 
transfer function (MTF) response at zero spatial frequency. As a 
starting point the voltage responsivity 9?/ of a single pixel can be 
expressed as 
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[V photon'^ 


(5-13) 


^ sec 

where Vj, is the voltage output of the device, and Ep is the uniform 
photon irradiance on the device [photon sec'^ cm"^] . If monochromatic 
irradiance is considered, the expression for the output voltage of 
a pixel can be rewritten as[29] 

Vc = Ep(X) Ax Ay t] e (CE) [V] (5-14) 

where AxAy is the resolution element area, is the integration 

time, r|, is the quantum efficiency, e is the electronic charge and 
CE is the conversion efficiency of the output preamplifier 
(V/electron) . The quantum efficiency is defined as the probability 
that a photoelectron is produced when a photon is incident on the 
detector. Substituting Vp in the responsivity expression Eq. (5-13) , 
the spectral array photon responsivity can be expressed as 

^ ^integ ^ ® [Vphoton'^cm^ssc] . (5-15) 

The responsivity, in Eq. (5-15) is valid for a pixel or 

for an array at zero spatial frequency. To include the resolution 
effects would require that the responsivity would be multiplied by 
the MTF of the array in two orthogonal directions, thus getting two 
expressions for spectral array photon responsivity in the x 
and y directions. 

The noise equivalent quantum flux (NEQF) is the counterpart of 
the NEP for a focal plane array. The NEQF is defined as 
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NEQF ■ 




(5-16) 


[photons sec'i cra'^] 


where NEQF has units of photon sec"'- cm"^ and noise is in V. 

The voltage noise of the device is caused by all the 
contributors shown in Fig. 5-1. A detailed discussion of the 
detector noise can be found in references[22][27][29]. Hence 




[V] 


(5-17) 


where An„^ is the sum of all the noises in quadrature in units of 
electrons (statistical independence). 

Substituting Eg.(5-15) and Eg.(5-17) in Eg.(5-16) one finds 


NEQF . 


Ax Ay h 


[photons 


(5-18) 


This is the photon irradiance flux on the array that would give a 
signal-to-noise ratio of 1, and is a measure of the minimum photon 
irradiance level that can be detected. The threshold signal-to-noise 
level will still be required for detection. The value of p will 
depend on the operating conditions. 
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Figure 5-1. Noise sources in a IR seeker using a CCD. 

In the case of background-limited infrared photodetection 
operation (BLIP) the dominant noise is shot or G-R noise associated 
with the photon irradiance. One can then write the NEQF as, 


NEQF . 


^2G^DE^Ax hy t 
Ax Ay r| 


(5-19) 


NEQF ■ G 


(5-20) 


where G is the photoconductive gain (G=l for photovoltaic operation, 
and 0<G<1 for photoconductor operation) , and t^^^es is the integration 
time. 

A photoconductor which exhibits performance at theoretical 
limits is said to operate in the BLIP mode. In this case the 
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background photon flux and the atmospheric transmission determine 
the detector noise. In this circumstance the flucuations in incident 
photon flux can be the dominant source of noise. For the case of an 
IR seeker where most of the field of view contains background and 
not signal photons, and where the background photon flux due to dome 
emission is significantly high, the term BLIP is truly appropriate. 
This is further enhanced by the fact that the IR sensor, in current 
systems, is cooled well below ambient temperature. 

B. SIGNAL-TO-NOISE RATIO DEGRADATION DUE TO DOME EMISSION 

In a high speed flight situation (^^,>1.2) the missile heats up 
due to aerodynamic heating. As was discussed in Chapter I, the 
temperature rise is particularly severe in the stagnation regions as 
happens with the IR dome. As the dome heats up, it will emit gray- 
body radiation over the entire electromagnetic spectrum. This will 
degrade the IR sensor performance in terms of detection and tracking 
ability. The effect of background-induced noise from the hot dome on 
the sensor performance has been characterized and explained by some 
authors, in particular Dr. Claude Klein. This section and the 
nonimaging range equation subsection summarizes some of his work. 

With the assumption that the seeker system has a wavelength- 
independent optical efficiency over the spectral bandpass, that the 
sensor subsystem is detector-noise limited, and that all detector 
elements operate in a BLIP mode, the peak-signal to rms-noise ratio 
obtained at the output terminals of the signal-processing unit can 
be expressed as follows[4]: 
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WEI 


(5-21) 


Here the symbol [W/cm^] designates the net effective irradiance 

originating from the target, 

^eff ■ f (5-22) 

where ip is the path transmittance and is the spectral irradiance 
originating from the target and measured at the collecting aperture. 
NEI is the noise-equivalent irradiance, 

WEI ■ WEi^yST^ (5-23) 


referred to the peak-response wavelength of the detector package, A^, 
and expressed in terms of the laboratory-measured NEI^ and the 
background photon-flux ratio 0/4),,. Under laboratory conditions there 
is no sunlight, and the system is in thermal equilibrium (the seeker 
and its surroundings are at the same temperature). Under background- 
limited conditions (BLIP) , the noise voltage is seen to be 
proportional to the square root of the flux of incident background 
photons(Poisson photon statistics), so that the signal-to-noise 
ratio (SNR) in the presence of hot IR dome radiation can be derived 
from the "dark-system" SNR simply in the form 


SNR . 




(5-24) 
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The ratio thus determines the performance degradation induced 
by environmental factors and must be evaluated with some care. Let 
one consider that the noise-generating photons originate from three 
distinct sources: the detector surroundings, the background scene, 
and the hot dome. In terms of spectral quantum irradiances, the 
total effective photon-flux density C> thus amounts to 

O - f (^x,B • *x,o • (5-25) 


where the terms in the integrand refer to the spectral photon flux 
contributions from the background, dome and detector surroundings 
respectively. The temperature distribution on the dome in practical 
designs does not vary significantly. So the dome can be considered 
to be at a uniform, averaged temperature . The spectral region of 
interest usually is narrow enough, implying that the emittance does 
not change significantly, so the photon flux can be expressed in the 
following compact manner[4], 


-1 


ri - e; CfV • ^ 

cTiy 


(5-26) 


Here F represents average dome emissivity, is the transmission of 
the optics(about 50% for state-of-the-art seekers), Q(T) represents 
the effective quantum exitance at absolute temperature T defined as 

Q. j dK . [photons/cm^-sec] . (5-27) 
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Under laboratory-type conditions, the effective background photon- 
flux density received by a radiation-shielded detector amounts 
to[26] 

. (NA)^ j (TJ dX , (5.28) 

Here NA is the numerical aperture of the optical system and Zf 
represents the (cold) filter transmittance,assuming ideal spatial 
shielding and spectral filtering. 

C. TARGET ACQUISITION RANGE DEGRADATION DUE TO DOME EMISSION 
1. Range Equation for Non-Imaging Systems 

A point target of spectral radiant intensity gives rise to 
a spectral irradiance 


(5-29) 


at the collecting aperture of the seeker system, where Zp represents 
the overall path transmittance, including the IR dome losses, and R 
is the target-collector range. Hence, and in accordance with Eqs.(5- 
20) and (5-21), we have 


j (X/\) J 


(5-30) 


where refers to a net effective apparent target radiant 

intensity defined as follows: 
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(5-31) 


The effective radiant intensity varies with slant range and IR dome 
temperature through the path transmittance Target acquisition at 
a given probability of detection and a fixed false-alarm rate 
requires a minimum signal-to-noise ratio (SNR)^^^. It follows that 
the detection or acquisition range is a solution of Eq.(5-29). If R 
now designates the acquisition range for the "cold" system obtained 
under laboratory-type conditions, the degradation induced by heating 
up of the IR dome can be characterized by[25] 


(5-32) 


(neglecting path transmittance 
from Eq.(5-31) 

(SNR) I 

so that using Eq.(5-21) , 

_ 5 . - 1 



Using the relationship 


(5-33) 


(5-34) 


since is a relatively weak function 

scenarios of practical interest. 


of both R and , for 


87 




FLUX DENSITY (3to5 pnoionsW• sec oi aeiector) 



Figure 5-2 Photosaturation characteristics of a lead telluride detector 
subjected to background illumination originated from a hot magnesium 
fluoride window[26]. 


The dark-system NEI relates to design parameters in the 
following manner[25]: 


(NEI)^ ~ 


(X^) 


(5-35) 


where D'(Ap) refers to the peak specific detectivity of the sensor 
package, under conditions of thermal equilibrium, at room 
temperature, and rig is an electric response factor. 

Figure 5-2 shows calculated and computed minimum detectable 
signal power relative to minimum detectable signal power in the 
absence of window radiation[26] . This example shows the strong 




)f the dome emission 


detection and tracking. 




2. Range Equation for Imaging Systems 

The NETD Eq.(5-10) is valid either for a BLIP or a non-BLIP 
detector. In this section the NETD for a BLIP detector irradiated by 
a hot window will be derived. As will be shown the range equation 
for an imaging system will be NETD dependent. 

The theoretical limit for D* for a photoconductor is 


D-(X) . —— 


(5-36) 


where is the quantum efficiency and is the background quantum 
flux density [photons/cm^'sec] incident on the detector. The severity 
of the IR seeker signal-to-noise-ratio limited application dictates 
that a cold filter, with appropriate spectral response, should be 
placed immediately preceding the detector. The Qg is given then by 



(T^) dX , 


(5-37) 


where it assumed that the detector is a Lambertian receiver with an 
effective collecting solid angle of n, and where is the ideal 
effective angle to which the detector cold shield reduces reception 
of background radiation, and r^o is the transmittance of the cold 
filter. At this point, it is convenient to define the several 
individual transmittances that affect the sensor performance. The 
total spectral optical transmittance can be decomposed in the 
following product 
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(5-38) 


■ '^co,\ ■ ■^ac.x ■ CV 
where t^c,a the spectral transmittance of the uncooled optics and 
is the average spectral transmittance of the dome at an 
averaged absolute temperature 

In this analysis it is assumed that the detector is a 
photoconductor having theoretical performance such that 


for X 

> K ■ 


Xp 


Substituting for Qg in the D" Eq. (5-36) gives 


(5-39) 


f ^co 2x 

A normalized detectivity parameter, D’‘ggj.p(A) , for unity quantum 
efficiency and for an effective cold shielded angle of n steradians 
is defined by[22] 

o'blip (>-> ■ T- - -TTr - 




Expressing D‘(A) in terms of D" (A) 
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(5-47) 




where Q^s = a-b / f ^. 

The perfect cold shield is defined as one which limits 
reception of background radiation to the cone defined by the optical 
system aperture ratio, F„ (f/n\mtiber) . Thus for an optical system with 
circular aperture of diameter and focal length f 


Qp . 


nn^ 

4f^ 


(5-48) 


and substituting in Eq.(5-43)gives[22] 


NETD . 


_ ^ _ 

y/oP' ^ 


(5-49) 


Applying the same definition to D‘ gives 

D- (\) . 2 Fn fri„ ■ (5-50) 

The Eq.(5-49) can be rewritten 


NETD . 


4 Fn 

ya|5 D-(X^) ^ 


(5-51) 


To account for the effect of the hot dome on the NETD the 
background quantum flux density from the dome emission should be 
added to the background quantum flux density from the scene 
background in the D’ Eq. (5-36) . With the cold shield, the background 
quantum flux density that irradiates the detector has two 
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contributions: 


2b ■ 


2d 


(5-52) 


where is the quantum flux density originated by a uniform 
blackbody background corresponding to the scene at temperature Tgg, 
Q„ is the portion of the background quantum flux density originated 
by the dome with a uniform mean tenperature . These quantities are 
defined as follows [25]: 


/ 


m / 


with radiometric constant as defined before and Cj = 

1.8837xlO^^Um^/sec'cm^, e is the mean dome emissivity with an 
averaged temperature, . Aj and A 2 are the cut-on and cut-off 

wavelengths of the filter. Substituting in Eq.(5-38) yields 


D-(\). 


X 

2nc 



(5-55) 


which can be used to express NETD from Eq.(5-50) in terms of photon 
flux components. 

Typically the optical system installed in the seeker is a 
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modified Cassegrain telescope. This will introduce some energy 
losses characterized by the circular aperture attenuation 
coefficient, K(cir), and a Cassegrain obscuration coefficient 
K(cas). If these losses are considered, the general NETD becomes 


/t [MbsiO-sb) ■'■Mb i(TB)]dA\ 


/cTp K(cas)K(cir) \ 


Tig - 


where to,- is the cold optics transmittance (under laboratory 
conditions) , to is the optics transmittance (during operation) , Mbb,x 
and Mg }, are the spectral radiant exitances of the background and of 
the dome respectively. Expanding the terms of the integrands of 
Eg.(5-56) another format to the general NETD equation can be found. 



(5-57) 

The last two expressions for the NETD take into consideration the 
influence of radiation by an aerodynamically heated dome on the 
sensor. The added "background"-induced noise will lower the 
effective performance of the system in terms of detection range and 
imaging capability. 

The signal-to-noise ratio at maximum detection range is 
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AT -TpCR; 

wire 


(5-58) 


where Tp(R) is the total atmospheric path transmission for the system 
spectral bandpass and R is the detection range. Here AT is the 
differential temperature of the target with respect to background. 
If AT is normalized with respect to the SNR, defining[30] 


Differential Tempen 


AT 


(5-59) 


and substituted in Eq.(5-58) then after rearranging one obtains 

■ —z — • 


To solve for range, the total transmittance of the optical path 
should be determined first. For a given set of meteorological and 
oceanographic conditions, the transmittance falls exponentially with 
range, as follows 


where a and 3 are constants to be determined. Substituting Eq.(5- 
61) in Eq.(5-60), gives 


95 




NETD 


(5-62) 


Solving for R yields the range equation for an imaging system 
In a - In(NETD) . In(Kl 

R ■ -p-^ . (5-63) 

in the form utilized by S. Smith[30]. 

D. DETECTOR NONUNIFORMITY NOISE 

This noise is associated with the variation in output from 
pixel to pixel, or the variations across the entire array. It can be 
caused by two effects:(a) fixed pattern noise and (b) responsivity 
variations. The influence of nonuniformity on FPA performance may 
impair the detector of achieving theoretical (BLIP) sensitivity. 

Nonuniformity limitations on the NETD in many systems are more 
damaging than temporal noise since a pattern-noise-limited image 
cannot be improved by frame-to-frame averaging, whereas the 
sensitivity of a temporal-noise-limited imager can be improved by 
such averaging. 

1. Fixed Pattern Noise 

Fixed pattern noise is a nonuniformity appearing across the 
array which does not change with irradiance[29]. Several effects 
contribute to fixed pattern noise. Thermally generated carriers are 
not uniform across the array. Local defects in the semiconductor 
crystal lattice cause a variation in dark-current generation across 
the array. 

A major contributor to fixed pattern noise is the errors in the 
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fabrication process. The errors introduced by inaccurate pixel areas 
caused by masking misalignments can lead to 10% nonuniformities. The 
fixed pattern noise is usually compensated by image processing. 

2. Responsivity Variation 

Responsivity variation is a nonuniformity which has a variable 
effect as a function of irradiance. The responsivity variation from 
pixel to pixel is determined by doping variations across the device. 
Compensation for responsivity is more difficult then for fixed 
pattern noise- requiring "two point" calibration as a minimum. 

3. Residual Nonuniformity 

After compensation the residual non-uniformity can be 
considered as a source of noise. This noise can be considered as an 
added percentage of photon noise, because it is considered to be 
correlated, in the following manner, 

(5-64) 

where n is the rms value of the noise in volt, g is the non¬ 
uniformity coefficient, Q^b and are the background photon fluxes 
on the detector from the scene and seeker dome as defined before by 
Eq.(5-53) and Eq. (5-54) respectively. Array non-uniformity is a 
problem area that has been discussed extensively in the literaature. 
For state-of-the art sensors the array non-uniformity coefficient 
varies between 0.01 and 0.0001(30]. 


E. COMPUTATION OF THE DETECTION RANGE 

Acquisition range for an IR seeker of an air-to-surface anti¬ 
ship supersonic missile with a 256x256 CCD staring sensor was 
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COMPUTATION OF THE DETECTION RANGE 


Acquisition range for an IR seeker of an air-to-surface anti¬ 
ship supersonic missile with a 256x256 CCD staring sensor was 
computed for a midlatitude scenario (LOWTRAN 6) , in winter and 
summer, for two background-to-target normalized differential 
temperatures (iCj^2.5 and JCj^=4.0) and for two missile flight profiles, 
level flight at 20 m and 3km altitude, and for a range of mean dome 
temperatures of 300, 600, 900, 1043 and 1200K, corresponding to 

missile flight speeds up to Mach 4. 

The detection or acquisition range was computed using Eq.(5-63) . 
To compute the NETD a MS-FORTRAN program developed, by S.T. Smith 
(Naval-Air Warfare Center, China Lake) was used. The program 
computes the NETD for different detector materials, integration 
times, dome materials, dome temperature, array nonuniformities, fill 
efficiency, background temperature, collector aperture and optical 
configurations. The NETD was computed for each set of parameters 
shown in Table 5-1, with LOWTRAN atmospheric transmission. 

The total transmittance of the optical path was determined 
first. The LOWTRAN6 computer code was used to predict atmospheric 
transmittance. Figure 5-3 shows the total computed transmittance 
for a slant path from 3km altitude to sea level and for a horizontal 
path at 20 meters altitude for several ranges in a Navy Maritime 
environment(see Appendix B) in the 3.5-4.2 pm and 8.0-10.5 pm, for 
the midlatitude summer and winter models. 
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Figure 5-3. LOWTRAN 6 computed atmospheric transmittance. * 3.5-4.2pm + 8.0-10.5 pm. 
Navy maritime aerosol model. (a)(b) Midlatitude Summer;visibility 44.2 km :Humidity 76.2%; 
(c)(d) Midlatitude Winter;visibility 43.9km:Humidity 77.1%. 

Using a least-squares fit, the coefficients a and of Eq.(5- 
61) were computed for the four different scenarios for the two 
spectral bands. The curve-fit parameters for atmospheric 

transmission and the computed NETD for the each set of sensor 
parameters (see Table 5-1) and dome temperatures of 300, 600, 900, 
1043 and 1200K, were used to compute the seeker detection range by 
applying Eq. (5-63) . The selected dome materials were the sapphire 
and the diamond for the 3.5-4.2vim band and 8.0-10.5vun spectral band 
respectively. These materials were matched to the spectral bands 
where they exhibit better transmission properties. These materials 
are considered to withstand the severe thermal environments of 
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advanced, highly supersonic missiles. The temperature of 1144K was 
found in Chapter IV to be the mean dome temperature of the 
conceptual missile at Mach 4 at sea level. 

The normalized differential temperatures, K^2.5 and K^4.Q were 
chosen for a SNR of 2.0. This implies that we are looking at 
temperature differences of 5. OK and 8. OK between the target and the 
background. 

F. RESULTS AND DISCUSSION 

Figure 5-4 through Fig.5-7 show the computed NETD for 
"residual" array nonuniformity coefficient after compensation of 
0.1, 0.01, 0.001 and 0.0001 respectively for the two selected dome 
materials. Fig.5-8 through Fig.5-23 show the computed detection 
ranges for the seeker. 

Before starting a detailed analysis of the detection range 
computation results it should be noticed that the values of the 
computed NETD depend on the area of the collector aperture(see 
Eq. (5-10)) . So all the discussion relates to the nose shape and 
dimensions of the hemispherical nose configuration of Chapter IV. 
The quantitative results obtained cannot be immediately 
extrapolated to other seekers. However the results obtained show the 
trends and influence of the several system parameters on the IR 
seeker performance. A LRASM (Long Range Anti-Ship Missile) would 
enable the installation of a larger primary optical aperture 
resulting in better detection range and good resolution for imaging. 

Analyzing the results it is clear that the array nonuniformity 
has strong effect on the seeker performance. The array nonuniformity 
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is particularly important at high supersonic speeds when the dome 
temperatures become particularly high. For the case of a Mach 4 
missile the seeker becomes blind (detection range falls to zero) for 
residual array nonuniformities of 10 percent. For the other cases of 
improved array nonuniformity the detection range of the system 
improves considerably. The sapphire dome systematically shows a 
better performance than the diamond dome in terms of detection 
range. 

In terms of environmental influence on the performance the 
winter condition with a slant path scenario gives the best results 
in opposition to the summer scenario for a sea-skimmer missile at 
an altitude of 20m. 

The results show clearly that the designer should determine 
carefully the effect of array nonuniformity and the dome temperature 
reached during flight by the seeker. 

The model used does not take into account target-to-background 
contrast reflectance. A better model should include this effect if 
a better target-to-background discrimination is desired. This effect 
is expected to be larger in the MWIR spectral band than the LWIR 
spectral band. 

For aim-point considerations a minimum of 100 pixels on target 
is desired (Johnson criteria). To calculate the aim-point tracking 
range one should take into account the geometrical characteristics 
of the target and the instantaneous field of view of the seeker 
optical system. 
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FPA.256x256, 30 pm pitch 

Aperture Diameter.2.25 cm 

FOV.4.46° 

IFOV.0.2° 

Diffraction limit blur.0.13mrad 

Obscuration coefficient.0.1 

Dome Material 

3.5 - 4.2pm band ( Sapphire ) 

8.0 - 10.5pm band ( ZnS and Diamond ) 


Detector Material 

3.5-4.2pm band .InSb 

8:0 - 10.5pm band.HgCdTe 


Array Nonuniformity 

0.001,0.01,0.10,1.00 

Dwell Time 

Staring.4.0 ms 

Dome shapes 

hemispherical 

elipsoidal 


Table 5-1. Seeker Parameters. Study Cases. 
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256x256 FPA 



Dome Temperature, [K] 
Figure 5-4 


256x256 FPA Integ. Time: 4msec 



Dome Temperature, [K] 






256x256 FPA Integ. Time; 4msec 



256x256 FPA Integ Time: 4 ms Kt=4.0 
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256x256 FPA 


iteg. Time; 



Dome Temperature, [K] 


Figure 5-9 






Kt=4.0 


256x256 FPA Integ. Time: 4 ms 



256x256 FPA Integ. Time; 4 ms Kt=4.0 
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256x256 FPA Integ. Time: 4 ms Kt=4.0 



256X256 FPA Integ. Time: 4 ms Kt=4.0 
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256x256 FPA Integ. Time: 4 ms Kt=4.0 



256x256 FPA Integ. Time: 4 n 



Dome *lemperature, [K.] 
Figure 5-15 
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256:256 FPA Integ. Time: 4 ras Kt=4.0 

Midlatitude Witifer 
honzontal patn,2u[n 
0.01 percent NonuniL 


400 600 800 1000 1200 


256x256 FPA Integ. Time: 4 ms Kt=4.0 


Midlatitude .Winter , 
slant path 3lcm-> sea level 
0.01 ^rcent Nonunif. 
Diamond 

■ - &0-10.5 microns 
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256x256 FPA Integ. Time: 4 ms Kt= 4.0 

80 
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le«:l 

0.001 percent Nonunit 


t>ome Temperature, [K] 
Figure 5-18 


256x256 FPA Integ Time: 4 ms Kt=4.0 



Dome Temperature, [K] 
Figure 5-19 
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Kt=4.0 


45 L -I 

I 0.00] percent Nonunit 


400 600 800 1000 1200 


0.001 percent Nonunit 
Diamond 
8.0-10.5 microns 


400 600 800 1000 1200 

Dome Temperature, [K] 

Figure 5-2 J 
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Dome: Sapphire ( £=1x10'^) 
Detector: InSb 

Spectral Band: 3.5-4.2 pm 

Dome: Diamond ( e=5xlCr^) 
Detector: MCT 

Spectral Band: 8.0-10.5pm 

Dome 

NEI 

NETD 

NEI 

NETD 

Temperature 

[Kl 

[W/cm^] 

[K] 

[W/cm^] 

[K] 

300 

3.293x10-'* 

2.236x10-' 

3.545x10'* 

1.026x10' 

600 

3.958x10'* 

2.689x10-' 

5.606x10'* 

1.623x10' 

900 

7.335x10-'* 

4.982x10' 

9.336x10'* 

2.702x10-' 

1043 

9.705x10'* 

6.592x10' 

1.140x10-" 

3.301x10-' 

1200 

1.270x10-"' 

8.628x10' 

1.379x10" 

3.993x10' 


Table 5-2. ComputedNEIandNETD for a 256x256FPA with 0.01% Nonuniformity. 



Dome: Sapphire ( e^lxlO'^) 
Detector: InSb 

Spectral Band: 3.5-4.2 pm 

Dome: Diamond ( £=5x1(7*) 
Detector: MCT 

Spectral Band: 8.0-10.5pm 

nnm. 

NEI 

NETD 

NEI 

NETD 

Temperature 

fK] 

fW/cm^] 

[K] 

[W/cm*] 

[K] 

300 

6.447Ixia'^ 

4.3790x1(7* 

3.2358x1(7'* 

9.3677x1(7* 

600 





900 

2.7617x1(7'^ 

1.8758x1(7' 

9.0180x1(7'* 

2.6107x1(7' 

1043 

4.4396x1(7'^ 

3.0155x1(7' 

1.1082x1(7'* 

3.2083x1(7' 

1200 

6.8278x1(7'* 

4.6375x1(7' 

1.3473x1(7'* 



Table 5-3. Computed NEI and NETD for a 256x256 FPA with 0.1 % Nonuniformity. 
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Dome: Sapphire ( e^lxlOt^) 
Detector: InSb 

Spectral Band: 3.5-4.2 pm 

Dome: Diamond ( e=5x1 (L^) 
Detector: MCT 

Spectral Band: 8.0-10.5pm 

Dome 

Temperature 

m 

NEI 

[W/cmf 

NETD 

[K] 

NEI 

[W/cm^J 

NETD 

[K] 

300 


3.789x10' 


9.358x10' 

600 


5.570x10' 


1.531 

900 


1.818 


2.609 

1043 


2.9574 


3.207 

1200 


4.5795 


3.870 


Table 5-4. ComputedNETD for a 256x256 FPA with LOO % Nonuniformity. 



Dome: Sapphire (e=lxl0'^) 
Detector: InSb 

Spectral Band: 3.5-4.2 pm 

Dome: Diamond ( €=5xICr^) 
Detector: MCT 

Spectral Band: 8.0-10.5pm 

Dome 

NEI 

NETD 

NEI 

NETD 

Temperature 

[K] 

[W/cmf 

[K] 

[W/crn^] 

[K] 

300 

3.2933x10" 

2.2368x10-" 

3.5448x10"' 

1.0262x10-" 

600 

3.958exl0" 

2.6886x10-" 

5.6063x10'“' 

1.6230x10-" 

900 

7.3345x10-" 

4.9817x10" 

9.3362x10'“' 

2.7028x10" 

1043 

9.7048x10'^ 

6.517x10" 

1.1401x10-'" 

3.3007x10" 

1200 

1.2703x10-'’ 

8.6283x10" 

1.3779x10'" 

3.9930x10-" 


Table 5-5. Computed NET and NETD for a 256x256 FPA with 0.0001 % Nonuniformity. 
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Dome: Sapphire (e=lxl(T^) 
Spectral Band: 3.5-4.2 pm 

Dome: Diamond ( e=5xl(T^) 
Spectral Band: 8.0-10.5pm 

Dome 



Temperature 

[K] 

HHHHH 

hmhhmiii 

300 

7.9172x10“ 

1.1486x10" 

600 

1.6141x10“ 

1.8792x10" 

900 

3.8033x10“ 

3.2031x10" 

1043 

6.1885x10“ 

3.9366x10" 

1200 1 

9.5829x10“ 

4.7860x10" 


Table 5-6. Computed Photon Incidence in the focal plane. 
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VI. COMPUTATION OF PHOTON NOISE FROM HOT DOME 


The characterization of fixed pattern-noise is an important 
aspect in evaluating the performance of a staring array. The hot IR 
dome will induce on the sensor a background photon noise 
distribution that can be considered like a fixed pattern-noise. In 
discussing the "hot dome problem" the question arises as to what 
dome temperature distribution should be used in the calculation of 
the dome emission effects on a staring sensor performance. Is it 
sufficient to consider the temperature of the dome to be uniform, 
and use a mean temperature, or will the background photon flux 
induced by dome emission impose a significant fixed pattern noise on 
the detector array? A useful simplification of the problem would be 
to consider the photon flux constant through the detector area and 
consider an averaged dome temperature. 

In this chapter these questions will be investigated by 
computing the radiative interchange between the hot dome and the 
detector array. 

A. RADIATION TRANSFER 

Given the temperature distribution on the nose, we are able to 
compute the IR emission due to the aerokinetic heating of the 
window. In our case we want to compute the self-emission generated 
by the hot window that will fall on the seeker's IR sensor. As 
assumptions, we will consider that the window has reached thermal 
equilibrium, and so the temperature of the dome is the recovery 
temperature, that the window is thermally thin, and that the dome 
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inner surface is not coated with any material. In order to compute 
the window's self-emission sensed by the detector, an optical system 
must be assumed. In this case a small telescope, composed simply of 
a single germanium lens, is adopted. This optical design may seem 
unrealistic (no field lens or light pipe are used), but in fact, it 
can serve as a benchmark to compare some of the aero-thermal effects 
that degrade the system image quality due to dome emission. 

To start the radiant transfer analysis, the radiant power 
reaching an element of area dA^ of the detector surface due to the 
source element dA„ on window surface S„ should be computed first. 
For simplicity one should assume, for the moment, that no lens is 
present between the self-emitting window and the detector. The line 
joining the element areas, of length r^,,, makes angles of &„ and 6^ 
with the respective normals to the surface elements, as shown in 
Fig. 6-1. The radiant power is d^W^, a second order differential 
because both the source and receptor are elemental areas[11][21], so 
that 


LdA dA cosd cosd. 

d ^ ^ ^ (6-1) 

where L is the radiance at dA„. 

The total radiant power at the entire second surface due to the 
entire first surface is, by integration. 

. - LcosB cos&^dA^dA 

" • (6-2) 
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Figure 6-1. Configuration for interchange between 
two infinitesimal elements. 

In adding power rather than amplitudes in this integration, we 
have tacitaly assumed that the radiation source emits incoherent 
radiation and that the radiation emitted from any point of the 
surface is diffusely distributed, obeying Lambert's law. If we 
discretize the window and the detector in finite area elements, and 
consider each element to be isothermal, the total irradiance into 
the detector area element originated by the window is 

J L.COsQ^ COsQ^ 

where the window is discretized in J finite area elements and the 


det( 


K f; 


The radii 


:tral 





region of interest of the window finite element j is 







(6-4) 


where c is the speed of light in vacuum, h is Planck's constant and 
k is the Boltzmann's constant, and i'(X,T) the emissivity of the 
infinitesimal element. 



Figure 6-2. Configuration for interchange between two finite 
surfaces. 


As was shown in the previous chapter, when the interest is in 
the limits imposed by dome heating on seeker performance, a more 
relevant quantity is the photon exitance Q (photons s"’^ cm'^) from the 
window. The spectral photon exitance, (photons s"^ cm'^pm"^ ) , is 
obtained by dividing Planck's law by hc/A, 
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(6-5) 


The integrated photon exitance,in a given wavelength band, is then 
obtained by integrating Eg.(6-5) between the cut-on and cut-off 
wavelengths, 


QfT) . 



(6-6) 


Numerical approximation of this integral over a narrow wavelength 
range yields 


QCr) 


2nc 


(6-7) 


where is the average wavelength within the band and A^. 

Then the photon fliox that irradiates a finite area element k of 
the detector is 


j E.(K,T)Q cosd cos6 A 


(6-8) 


where Qj 


black body exitance of th 


ilement j. 





COMPUTATION METHOD 


Now let's consider that between the window and the detector 
there exists a thick converging lens. This system is not in thermal 
equilibrium. In fact for simplification consider that the detector 
and lens are cooled to zero absolute. To compute the irradiance on 
the detector due to the thermal emission of the heated window, a ray 
tracing must be done. In this method, the paraxial system matrix is 
used to compute the image surface for the lens conjugate to the hot 
irdome (infrared dome) as object surface. This image surface is then 
treated as an equivalent input surface viewed directly from the 
detector plane. The conjugate surface is found by mapping several 
points of the window by means of a matrix transformation. This 
virtual surface is then divided into diffusely emitting finite 
elements. Radiative transfer calculations between each detector 
element and the finite elements of the image surface are carried 
out, taking into consideration the respective shape factors. The 
contributions of all finite elements are integrated, enabling the 
computation of the photon flux that irradiates each photo-sensitive 
detector element. This procedure is applied to all the elements of 
the detector array. An array fill factor of 1.0 is assumed. 

1. System Ray-Transfer Matrix 

If we consider the temperature field of the window to be 
axisymmetric, resulting from zero incidence of the missile, and the 
optical system to be looking in the missile longitudinal axis 
direction, the problem of computing the photon flux that falls into 
the detector becomes a 2-dimensional calculation with cylindrical 
symmetry. Considering the optical system consisting of a single 



lens, a matrix can be constructed representing the action of the 
thick lens on the radiation emitted from the window, as follows: 


M 


S3 *3 S3 Sj Sj 


(6-9) 


where M the composite transfer matrix is formed by: translation 
matrix, 3 1 , representing the optical path from the input plane to 
the front face of the lens; refraction matrix, , representing the 
refraction process at the lens front surface; 3 2 < translation 
matrix through the lens medium; 9?2 , describing the refraction at 

the lens back surface; and finally 3 3 , describing the final 
optical path from the lens back surface to the output plane, 
conjugate to the input plane. 

The system matrix has the following format. 


M - 


(6-10) 


If the output plane corresponds to the conjugate to the input 
plane, the matrix element must be identically zero [21] . Then 
is the transverse linear magnification. 

2 . Calculation of the System Ray-Transfer Matrix 
In this section the system matrix will be calculated by 
multiplying consecutively the several translation and refraction 
component matrices for paraxial rays. 
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(6-1 la) 


(6-1 lb) 

(6-1Ic) 

(6-1 Id) 

(6-1le) 


where is the distance from the input surface to the lens, and 
i ?2 are the radii of curvature of the front and back surfaces of the 
lens, n and the indexes of refraction of the air and of the lens 
respectively, t represents the lens thickness, and finally, x is the 
distance from lens to the output plane. Carrying out the matrix 
multiplication as in Eq.(6-9) yields the following matrix elements 
and of the system matrix M: 
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3 . Computation of the Footprint at the Dome 

The computation of the footprint at the dome of the system is 
necessary in order to determine which portion of the dome will 
irradiate the detector, and so contribute to the background noise. 
This can be done by tracing back the limiting optical rays from the 
detector that pass through the lens aperture. Again considering the 
same system geometry, the problem becomes 2-dimensional. We have 



(6-15) 


where is the angle subtended by the ray from the detector passing 
through the lens margin, as shown in Fig. 6-4, and the 
composite matrix M‘ defined as follows: 

M‘ . S‘2 St) S'l (6-16) 

where , Sj' =32 - 3^2' = ^2 < a-nd the translation matrix 
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detector into the !< 


back surface is defined as folic 


(6-17) 


The matrix multiplication as written in Eq.(6-14) leads to the 
following matrix M' elements 




(6-18a) 





(6-18c) 





(6-18d) 


Substituting the following values into the respective expressions. 


n =1.000 
nJlOum) = 4.003174 
Ri = -50.000cm 
Rj = 60.000cm 
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t = 0.300 cm 
tj = 9.09099 cm 
Rap = 1.500 cm 






yields the following matrix: 


0.996249 9.131830 

-5.683183e-2 -1.483700 


with given by 


Using Eq.(6-13), 


y [o.996249 9.13183ol 0 

sinaj ■ [ .5.683183e-2 -1.483700 J sin(9.37°) 


yields 


(6-19) 


(6-20) 


(6-21) 


y 1.4S7 
sinaj ’ [- .2416] ' 


(6-22) 


which implies that a = 13.98° Furthermore, it follows that 


(z-rj^ . . ri (6-23) 

where r„ = 3.3166 cm (dome radius of curvature) . This defines a 
circular pupil on the dome of radius r = 1.9 cm (z=0.595 cm), 
computed by the following approximate formula: 

r . tga(t, - z) . . (6-24) 
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The quantities 


tj, 


and Rj, 


defined 


Fig. 6-4. 





Figure 6-4. Ray tracing variable nomenclature. 


C. RESULTS AND DISCUSSION 

The characterization of fixed pattern noise is an important 
aspect in evaluating the performance of a staring array. Fixed 
pattern noise is a generic term used to describe the variations in 
gain and offset levels observed over the FPA detector elements. 

The motivation of our study is the characterization of 
fixedpattern noise induced by dome emission. 

Computations were performed for the case of a hemispherical 
dome corresponding to the case I nose configuration (see page 49) of 
a missile at Mach 4 at sea level with zero incidence angle. Two 
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different computations were performed, the first considering that 
the dome temperature field was that obtained by solving the thin- 
layer Navier-Stokes equations and considering an adiabatic wall as 
one of the boundary conditions; in the other case the dome 
temperature was considered to be the mean value of the first case. 


DOME TEMPERATURE 



Figure 6-5.Thin-Layer-Navier-Stokes computed 
dome temperature. 

The optical system was considered to be looking straight ahead 
and just the portion of the window corresponding to the system's 
window pupil was considered. The collecting lens was designed using 
the ZEMAX XE ray tracing program. The lens data are given in 
Appendix E. The dome material chosen was CVD Diamond. The waveband 
considered was from 8.0-10.5pm. Fig. 6-5 shows the dome temperature 




versus off-axis angle. For the case considered the temperature field 
has cylindrical symmetry. Fig. 6-6 shows the photon irradiance on a 
64x64 FPA(Focal Plane Array) due to this temperature distribution. 
Fig. 6-7 shows the photon irradiance on a 64x64 FPA based on the 
mean dome temperature of the same temperature distribution. The mean 
dome temperature considered is 1024 K. 

For the computations cell separation of 30pm and a fill factor 
of 1.0 are considered. A fill factor of 1.0 although unrealistic 
serve the purpose of showing the variation of the photons flux 
through the array. 

The usual means of characterizing fixed-pattern noise is by 
specifying the variance, a, of the pixel levels, 


E E CiK-y,; - M]^ 

number of pixels 


(6-25) 


The variance computation is often written(schematically) as: 

. (I ■ M)^ (6-26) 


the array data being denoted by I(x,y) and the mean by M. A 
relatively larger variance is seen for arrays exhibiting an 
appreciable amount of fixed-pattern noise. However, a specification 
of variance alone gives no indication as to how the non-uniformities 
are distributed over the array, essentially assuming that there is 
an equal contribution from the fixed-pattern noise at all spatial 
frequencies. 

Figure 6-6 shows the percentage variation of the quantum 
irradiance on a 64x64 focal plane array due to a dome with a non- 
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uniform temperature field distribution. The values were computed in 
relation to the maximum quantum irradiance at the center of the 
array. The maximum quantum irradiance is 9.924x10’-^ photons/sec per 
pixel at the center of the array. The minimiim quantum irradiance is 
9.839x10’-^ photons/sec per pixel at the corners of the array. The 
computed mean quantvim irradiance is 6.2829 photons/sec per pixel. 
For this case the maximum computed irradiance variation is -0.8% at 
the corners of the array. 

Figure 6-7 shows the percentage variation of the quantinn 
irradiance on the focal plane array due to a dome with a uniform 
temperature field distribution. The values were computed in relation 
to the maximum quantum irradiance at the center of the array. The 
minimum quantum irradiance is xlO^^ photons/sec per pixel at the 
corners of the array. The computed mean quantum irradiance is 
6.2829 photons/sec per pixel. For this case the maximum computed 
irradiance variation is -0.8% at the corners of the array. 

This is five times more than the variance found for a uniform 
temperature dome at the average computed temperature. (Fig.6-7). 
The computed variance is 1.26x10^^photons^ per pixel for the case of 
a dome with uniform temperature. 

Fig. 6-8 shows a detectivity contour plot for the case of 
uniform dome temperature. For that purpose the following equation 
[29] was used. 



where Qp® is the background photon flux per pixel element area. The 
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Figure 6-8. Detectivity contour on a 64x64 FPA induced 
by the fixed pattern of the dome background irradiance 
(non-uniform dome temperature) 
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VII. CONCLUSIONS 


As expected, flight at high supersonic speeds has a dramatic 
effect on target detection range degradation of IR guided missiles 
due to the aerokinetic heating of its IR window dome. The 
temperature increase of IR dome increases the background quantum 
flux density that irradiates the IR sensor leading to an decrease of 
the signal-to-noise ratio and the correspondent decrease of the 
NETD of the detection system. An equation was derived that shows 
this effect. For Mach number speeds of 2 and 4 the detection range 
degradation for the cases considered reached values around 37% and 
97% respectively for array non-uniformity of 1%. This means that at 
high supersonic Mach numbers the seeker HR ("Imaging Infra-Red") 
system becomes "blind" due to dome emission. For lower array non¬ 
uniformities, i.e. 0.001 percent the detection range degradation was 
smaller in the range 12% and 22% for Mach 2 and 4 respectively. 

Since the integration times of the array are much smaller than 
the missile guidance system time constant (over 100 times) the 
photon irradiance due to dome emission on the FPA originates a 
fixed pattern noise on the detector responsible for a detectivity 
non-uniformity of the array pixels. As shown in Chapter 5, the IR 
seeker detection range has a strong dependance of the detection 
array non-uniformity. 
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APPENDIX A. CALCULATION OF DOME THICKNESS 


Any viable missile dome structure must be able to withstand the 
mechanical stresses induced by aerodynamic pressure as well as 
thermal stresses induced by aerodynamic heating. In a thermally thin 
regime, the TSR (Temperature Shock Resistance) performance will be 
maximized by making the dome as thin as possible. The TSR indicates 
the range of Mach ntunber variation that an IR dome can withstand 
without a good chance of catastrophic failure[10]. This temperature 
corresponds to the maximum difference in temperature between the 
stagnation temperature and the temperature of the inner wall of the 
dome that the dome can withstand. This leads to an ''optimum" dome 
thickness equal to the minimum dome thickness required to 

withstand the aerodynamic loads. Since the emittance of the dome 
increases with its thickness the minimum dome thickness will 
minimize dome self-emission effects as well. An accurate prediction 
of the dome emission effects necessitates an accurate computation of 
the dome temperature and of its emittance, and by consequence of the 
dome thickness. 

1. Pressure Induced Stress 

In a supersonic flight environment, altitude and velocity 
determine the pressure load. We should follow the procedure outlined 
in section II.C in order to calculate the minimum dome thickness to 
withstand the aerodynamic load, 

fPst ,1 

Ap ■ -p_ - p. - 1 (A-1) 
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where 



Consider a medium altitude flight profile. From II.S. Standard 
Atmosphere Tables one has for an altitude of 3 km: 

Y=1.4 

p..,=10.2 psi . 

Evaluating Eq.(A-2) gives a pressure jump along the stagnation 
streamline of 21.07. Substituting values in Eq.(A-l), gives Ap = 
204.7 psi. 

Let us consider that a "floating attachment" scheme, like a 
polyamide adhesive, is used to decouple the dome from the missile 
body. The maximijm tensile stresses failure criterion leads to the 
following expression 


a/SF 

Ap 



(A-3) 


for the minimum dome thickness, L^„ , where refers to the nominal 
flexure strength of the material, and SF is an appropriate safety 
factor that yields the minimum thickness required to prevent 
catastrophic failure of the dome. Since currently available IR 
transmitting dome materials are brittle ceramics, they may exhibit 
a wide range of fracture probabilities under apparently identical 
loadings, which suggest safety factors of at least four (4) if Of is 
the characteristic flexural strength derived from a Weibull 
statistical analysis. The dome thickness required for withstanding 
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aerodynamic loads Zip with good probability of survival can be 
obtained from the following equation: 



R 



(A-4) 


For a truncated hemispherical dome with R = 3.3166 cm, and using 
CVD diamond as dome material. 


75 (3.32xl0-^m)\ 


l.AlxlO^Pa 
1.5x10^Pa 


5. 4xl0'^m 


(A-5) 


2 . Thermally Induced Stress 

As shown in Chapter 1 the maximum tensile stress should be 
given by an expression such as [10] 


(A-6) 


where a = thermal expansion coefficient, 

TSF = Thermal Stress Factor, 

E' = biaxaial elestic modulus, 

AT = temperature differential across the dome wall. 


The stress induced by the sudden aerodynamic heating of the dome 
depends on the ability of the dome to conduct the heat flow and 
ultimately upon the heat-flow pattern, which is controlled by the 
Blot number. 
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Bi ■ 


(A-7) 


hL 

k 


where h = heat transfer coefficient at the outer surface, 

L = dome thickness, 

k = theirmal conductivity of dome material. 

The "optimum" dome thickness, L = is determined by 

aerostructural requirements and should be as thin as possible to 
enhance the thermal shock resistance and reduce self-emisssion. In 
assessing the "ultimate" TSR of an IR dome, the Biot number should 
be expressed in the following way 
where h ’ refers to a 1-cm radius dome. 

At Mach=4,■altitude=3km 


0.908y/R 


(A-8) 


Bi- . 


(0.908) (3.31 66cm; 
(20Wcm-^K-^) (150)^^^ 


1x10* 


(A-9) 


where R is in centimeters, k in watts per centimeter per kelvin, and 
Of is in megapascals. Considering the hemispherical dome window case, 
R = 3.3166 cm and CVD Diamond as window material, we conclude that 
the response of the diamond to transient heating will be in the 
thermally thin mode (Bi <1) . Indeed this material exhibits extremely 
low Biot numbers(Bi<0.001) that minimize the temperature gradients 
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across the thickness. 

The thermal shock temperature (TST) of the minimum thick 
temperature missile domes with Bi’ that represents the Biot number 
for the most severe environment based on projected tactical 
requirements. 

For thermally thin domes the peak temperature differential 
between the outer and inner wall is given by 

(iST)^ - ■ T.JxBi , Bi<l . (A-10) 

In the thermally thin regime the temperature of the back face 
may get hot long before the front reaches the recovery temperature. 

Heat Transfer Coefficient at Stagnation 



Figure A-1. Heat transfer Coefficient at Stagnation. 


141 







APPENDIX B. PROGRAM TO EXTRACT FLOW VARIABLES 

The following program written in FORTRAN was used to extract 
from the g.save file (PLOT3D format) the values of the flow 
variables of Table 4-2. It uses the grid.in file (PL0T3D format) and 
the g.save file as input files. By a simple alteration of the source 
code it is possible to read from the g.save computed file flow field 
values for a specific range of Kstreamwise direction), 
J(circumferencial direction) and K(radial direction). 
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C Extract x,y,z vs. Temperature, Normalized Temperature, 

C Coefficient of Pressure from grid and Q files in PLOT3D format, 
C 

C XYPLOT and MATLAB formats. 

C 

PARAMETER (MI=300,MJ=100, MK=50) 

DIMENSION XYZ (MI,MJ,MK, 3) ,Q(MI,MJ,MK, 5) 

DIMENSION X(MI) , Y (MI) , Z (MI) ,T(MI) 

CHARACTER*80 FILE 
C 

PARAMETER (GAMMA=1.4) 


C Read grid file. 

C 

WRITE(*, 1) 

1 FORMAT(' Enter PLOT3D grid filename: ',?) 

READ(*,2) FILE 

2 FORMAT(A) 

OPEN (UNIT=1, FILE=FILE, STATUS='OLD' , FORM=' FORMATTED' ) 

READ (I,*) NI,NJ,NK 

IF (NI.GT.MI .OR. NJ.GT.MJ .OR. NK.GT.MK) THEN 
WRITE(*, 3) 

3 FORMATC Grid is too big, exceeds internal dimensions.') 
GOTO 100 

ENDIF 

READ (1, *) ( (( (XYZ (I, J, K,NX) , I = 1,NI) , J=1,NJ) , K=l, NK) , NX = 1, 3) 

CLOSE(UNIT=1) 

C 

C Read Q file. 

C 

WRITE(*, 4) 

4 FORMAT(' Enter PLOT3D Q filename: ',$) 

READ(*,2) FILE 

OPEN (UNIT=1, FILE=FILE, STATUS^' OLD' , FORM=' FORMATTED' ) 
READd,*) NID,NJD,NKD 

IF (NID.NE.NI .OR. NJD.NE.NJ .OR. NKD.NE.NK) THEN 
WRITE(*, 5) 

5 FORMAT(' Q file dimensions don''t match grid.') 

GOTO 100 

ENDIF 

READd,*) FSMACH,ALPHA, RE, TIME 

READ (1, *) ( ( ((Qd, J, K,NX) , 1=1, NI) , J=1, NJ) ,K=1, NK) , NX=1,5) 

CLOSE(UNITED 
C 

C Get range of (i,j,lc). 

C 

C WRITE(*,11) 

C 11 FORMAT(' Enter IS, IE, JS, JE, KS,KE; ',$) 

C READ(*,*) IS, IE, JS, JE,KS, KE 
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APPENDIX C. OVERFLOW INPUT/OUTPUT FILE SPECIFICATIONS 


The following list describes the input/output files and the data 
they contain: 



Unit 

I/O 

JOescription 

grid-in 

I 

input 

grid file (PLOT3D format) 

fort.2 

2 

output 

input from PEGASUS 

q.restait 

3 

input 

0 restart file (PLOT3D format) 

q.save 

4 

input 

0 save file (PLOT3D format) 

q.bomb 

4 

input 

0 bomb file (PLOT3D format) 

resid.out 

7 

output 

flow solver residual history 

fomo.out 

8 

output 

force and moment coefficient history 

rpmin.out 

9 

output 

minimum density/pressure/gamma history U 

turb.out 

11 

output 

1 

turbulence model residual history 

species.out 

12 

output 

species continuity residual history 

interp.n 

10 

both 

PEGASUS interpolation and boundary point indices 

param.n 

10 

both 

intermediate storage of run parameters 1 

flowvar.n 

10 

both 

intermediate storage of flow variables 

gridmet.n 

10 

both 

intermediate storage of grid and metrices 


Table C-J. OVERFLOW Jnput/Outpulfiles. 
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APPENDIX D. LOWTRAN ATMOSPHERIC PROPAGATION MODEL 

The Atmospheric Transmittance/Radiance Computer Code LOWTRAN 6 
calculates both the atmospheric transmittance and radiance 
throughout the infrared spectrum(0.25 to 28.5 pm) [28] . It uses a 
single parameter band model for molecular absorption, and includes 
the effects of continuum absorption, molecular scattering, and 
aerosol absorption. Refraction and earth curvature are included in 
the calculation for slant paths[28]. 

Five seasonal models and the 1962 U.S. Standard Atmosphere are 
provided as inputs to the LOWTRAN program using typical altitude, 
pressure, temperature, water vapor density and ozone density 
profiles for each. The five other atmospheric profiles consist of 
(1) Tropical (15° N) , (2) Midlatitude Summer (45° N-July) , (3) 
Midlatitude Winter (45° N-January), (4) Subarctic Summer (60° N) , and 
(5) Subarctic Winter (60° N-January). Additionaly the user may input 
radiosonde data if desired to make calculations for a specific 
profile. 

The LOWTRAN code includes three boundary layer aerosol types: 
Rural, Urban, and Maritime. Due to the nature of the scenario, only 
the latter was used to simulate an engagement of an anti-ship 
missile in open waters. These aerosols are largely sea-salt 
particles in the lower boundary layer caused by evaporation of sea- 
spray droplets. Together with a background aerosol of more or less 
pronounced continental characteristics they form a fairly uniform 
maritime aerosol which is representative of the lower 2 to 3 
kilometers of the atmosphere over the oceans. 
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These aerosols should not be confused with the heavy direct 
sea-spray aerosols found in the lower 10 to 20 meters above the 
ocean surface which are strongly dependent on recent wind and sea 
surface turbulence. The extinction and absorption coefficients are 
calculated as a function of the relative humidity and based heavily 
upon the refractive index for the given altitude and weather 
conditions. 

The ability to calculate slant-path transmittance is 
incorporated in the LOWTRAN VI utilizing a data base of seasonal and 
latitude dependencies of vertical aerosol distributions. This allows 
transmittance calculations of medium and high altitude missile 
attack profiles. 

The LOWTRAN code calculates transmittance as a function of 
"absorber density" for the path, the pressure, and temperature, for 
the particular wavelength band chosen. It utilizes both empirical 
laboratory data and available molecular line constant data in 
performing its calculations. 

The aerosol population found over the world's ocean is 
significantly different in coitposition and distribution from that of 
a continental origin. These aerosols are largely derived from the 
sea. They are produced by the evaporation of sea spray and from jet 
and film droplets. Jet droplets are ejected into the air by the 
bursting of small air bubbles at the sea surface. The bursting of 
the bubble film leaves behind many smaller film droplets that may 
also be diffused into the air. These mechanisms are wind dependent 
and require Whitewater phenomenon to produce aerosol. 

Once the aerosol droplets are airborne, they undergo additional 
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nd mixing processes. The marine boundary layer is usually 
a temperature inversion and, within the boundary layer, 
sr marine aerosol together with any background aerosol form 
uniform aerosol spatial distribution. Once introduced into 
phere the lifetime of an aerosol particle is dependent on 
of the particular aerosol particle. Those with very small 
a a very long residence time in the boundary layer if there 
shout processes taking place. On the other hand, those with 
e sizes have a short residence time and do not contribute 
tationary long-term aerosol population. The reader is 
ed to Reference[28] for more details. 
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Figure E-2. Seeker collecting lens geometric 



Figure E-3. Seeker collecting lens pc 











APPENDIX F. "POLE-FIT" EMISSIVITY MODEL 


The "classical pole-fit" model" represents the complex 
permittivity, e, in the form of the lattice vibrations[24]: 


_ he (T)v](T) 

e(v,T) ■ e_ . 2^ —- i - , 

j Vj(T) - . iy.(v,T)v 


(F-1) 


where Ae^ , are the mode strength, line width, and long wave¬ 
length transverse optical frequency, respectively. This model 
adequately describe the reflectance, R, in the infrared region. 
Using Eq.(F-l), the static dielectric constant, e^{T) = e{0,T ), 
becomes 

eJT) . ejT) . J^Ae.(T) . 

The frequency dependence of y^ represents a cutoff beyond the 
one-phonon region caused by anharmonicities of the potential. The 
cutoff is empirically represented as follows 

Yj ■ Yj (T) ■ 

-exp |-oi I 

where a and v are arbritary parameters. Michael Thomas team at APL 
has been most successful in this effort by choosing a = 4 and T as 
1.1 times the highest infrared allowed transverse optical optical 
mode frequency. 
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As Eq.(F-1) suggests, the reflectance is a function not only of 
frequency but also of temperature. The temperature dependence of 
the transverse mode frequencies, , and the mode strengths, Ae^, 

can be represented by linear equations of the form 

V„, (T) . (TJ - a^(T- (F-4) 

and 

Ae^. (T) . Ae. (TJ . b.(T - TJ (p. 5 ) 

where is a reference temperature and and bj are mode-and 

material-specific constants coefficients. The dependence of the 
linewidth on temperature can be described by a quadratic equation in 
T as 

(T) . (TJ . C.(T - TJ . d^(T - , (F.6) 

'^TOj '^lOj 

where is again a reference temperature and Cj and dj are mode 
and material-specific constant coefficients. The high-frequency 
permittivity, e„ , is also temperature dependent and is represented 
by 

e_ (T) . e_ <TJ . e^(T - TJ . (F-7) 

ej=de^/T=2n^(dii^/dT); thus, in most cases, is determined from 

dn/dT measurements at the helium-neon laser frequency. 
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APEENDIX G. RESULTS FROM ELECTROMAGNETIC THEORY 


For optically smooth surfaces that are physically and chemically 
clean, electromagnetic theory with the aid of Kirchoff's law 
provides' a prediction of the magnitude and distribution of the 
directional emissivity[12]. The radiation emission properties of the 
aforementioned ideal surfaces are sometimes distinguished by using 
the term emissivity instead of emittance. 

Radiant energy either incident on or leaving from a surface 
can be resolved into two components of polarization. One of these 
components is polarized parallel to the plane of incidence and the 
other is polarized perpendicular to the plane of incidence. The 
electromagnetic theory predicts the specular reflectance (f 
corresponding to each component of polarization, defined as the 
ratio of the monochromatic intensities of that component in the 
reflected and incident beams. 

The expressions for p/ (perpendicular polarized conponent) and 
P||° (parallel polarized component) constitute the well-known Fresnel 
law of reflection assuming that the medium adjacent to the 
reflecting surface has a refractive index of unity and an extinction 
coefficient of zero (vacuum and gases). These expressions are 
derived for monochromatic radiation. However, the subscript A is 
omitted in order to avoid complicating an already cumbersome 
notation. The resulting equations are 






The angle 6 is the specular-ray direction corresponding to an 
incident beam whose angle of inclination is W = 6 . The quantities 
a and b are related to the angle of reflection 9 and to the 
refractive index n and the extinction coefficient k of the material 
on which the radiation is incident. 
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If the incident radiation is unpolarized, then the two components 
of polarization in the incident beam are of equal intensity. In this 
case it is readily verified that the monochromatic specular 
reflectance is the average of and pf, that is, 

p=r8; - -^[pf(e) . p^ue) ] ■ (g-4) 

The monochromatic directional absorptivity, a{d) , can be determined 
using Kirchoff's law, 

a ce; . 1 - p * re; - I - -^[pf(Q) ■ p,^(Q)j (g-5) 

Due to the fact that Kirchoff's law applies directionally and 
monochromatically for each component of polarization, one can write 
the components 


sjB) . ajd) . 1 

- p! fe; 

(G-6b) 

E, re; . re; - i 

- pfre; 

(G-6a) 

Furthermore, the directional emissivity of 

the mixed radiation is 

the average of Sj. and Sj , 

that is. 


ere; - 

• e, re;; 

(G-7) 


Although the subscript X is ommited for notational simplicity, the 
foregoing equations, as well as those that follow, are to be 
regarded as applying monochromatically. 

The predictions of electromagnetic theory for p/ and pj® can be 
employed to provide an expression for the directional emissivity. 
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Upon substituting equations (G-1) into equations (G-6) e (G-7), one 
obtains 



cos^B (a ' 


(G8) 


where 


E^ce; . 


4a cosQ 

■ 2a sinQ . cos^B 


(G-9) 


The quantities a and b are related to the optical constants and to 
the angle of emission d by equations (G-2). If n and k are known, 
Eq.(G-8) provides a prediction of the directional distribution of 
the monochromatic emissivity of an optically smooth surface. 

The hemispherical emissivity e can be determined from the 
directional emittance by integrating over the hemisphere. By 
introducing the definitions of s and £{&), the emissivity can be 
computed using the following relation 

^ ■ 4 / ^ ® ■ (G-10) 

Eq.(G-8), can be then be employed to provide a prediction for e. In 
the case of an electric nonconductor (k = 0), the execution of the 
integration yields 
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1 [1. in f 

2 3 3n 2(n^ . 1)^ 

(n .1) (n^ .2n -1) 4n^ (n^ .1) ^ 

. 1)^ (n -1) (n^ . 1)^ (n . 1)^ 


(G-11) 


in which e„ is the emissivity in the direction normal to the surface 


(G-12) 


The ratio of the hemispherical to the normal emissivity ratio 
has practical significance inasmuch as it is the normal emissivity 
that is more often determined by experiment, whereas it is the 
hemispherical emissivity that is usually needed in the computation 
of the radiant interchange. The ratio e/e„ is presented in Fig. G-2 
as a function of the refractive index for parametric values of the 
extinction coefficient = k/n. Within the range of n and ;c, values 
that correspond to real materials, the £:/s„ ratio is greater than 
one for metals and less than one for dielectrics. For the latter 
materials, f is only slightly less then 
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Figure G-2. Ratio of the hemispherical to the normal 
emissivity as predicted by electromagnetic theory. 



APEENDIX H. RADIATION TRANSFER SOURCE CODE 

The program listed below was written in MATLAB and computes the 
background photon flux on an nxn FPA originated by a hemispherical 
dome. The method used follows the procedure outlined in the section 
A and B of this chapter. The program assumes that the dome 
temperature has cylindrical symmetry and that the optical system is 
looking straight ahead. 

1. Input Parameters 

The temperature distribution, in kelvins, is entered in the 
vector tempdata. The first element is the dome stagnation 
temperature.The other temperature information is chosen in points 
downstream in order to ensure that a reasonable coverage of the dome 
is achieved. The coord has the radial coordinates (measured from the 
dome axis in cm) of the points corresponding to the data of vector 
tempdata. 

Following is a list of the input parameters: 

-> NW - radial discretization number of elements of the 
dome ; 

-> MW - number of sectors used to discretize the dome; 

-> Rpupil - radius of the pupil measured from the dome 
axis(m); 

-> R1 - radius of curvature of lens front surface(m); 

-> R2 - radius of curvature of lens back surface(m) ; 

-> Rdome - radius of curvature of the dome(m); 

-> t - lens thickness; 

-> tl - lens separation from the dome; 


163 






> t3 - backsurface focal length; 

> lambdamin - cut-on wavelength(microns); 

> lambdamax - cut-off wavelength(microns); 

> ND - sqrt(number of array pixels); 

> cell separation(microns); 

> em - dome emissivity; 

> n - index of refraction of the media suurronding the 
lens ; 

> nl - index of refraction of the lens; 

> toptics - cold optics optical transmission. 
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